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1.  INTRODUCTION 


This  report  presents  the  results  of  a  detailed  analysis  of  regional  wave  attenuation  along  con¬ 
tinental  paths  to  the  NORESS  array  in  Norway.  Our  data  consist  of  stable,  array-averaged 
spectra  from  190  regional  events  recorded  by  the  high-quality  NORESS  digital  instruments. 
The  logarithm  of  the  amplitude  spectra  from  events  of  varying  magnitude  and  epicentral  dis¬ 
tance  are  inverted  simultaneously  for  both  the  source  moment  and  the  apparent  attenuation.  The 
result  is  a  simple  parameterization  of  the  observed  amplitude  spectra  that  can  be  used  to 
address  a  number  of  seismological  problems  related  to  wave  propagation  in  the  region  and  to 
the  treaty  monitoring  capabilities  of  small  regional  networks. 

1.1  Objectives 

The  overall  objective  is  to  systematically  characterize  the  spectra  of  regional  phases  recorded  at 
NORESS.  The  emphasis  is  on  the  separation  of  source  and  path  contributions  to  observed  sig¬ 
nal  amplitudes.  To  accomplish  this,  we  have  developed  a  least-squares,  generalized  inversion 
that  simultaneously  estimates  source  strength  and  Q(f).  The  data  are  parameterized  by  an  (D2 
source  spectrum  with  cube-root  corner  frequency  scaling  and  an  assumed  geometric  spreading 
function.  The  method  utilizes  both  the  spectral  and  spatial  decay  of  observed  signal  amplitudes 
to  separate  source  and  path  contributions.  The  result  is  accurate  attenuation  curves  for  the 
regional  phases  analyzed.  Important  applications  of  these  curves  include: 

•  Simulation  for  hypothetical  sources.  Once  the  path  contribution  has  been  isolated, 
frequency-dependent  signal  amplitudes  resulting  from  a  known,  theoretical  source 
can  be  predicted.  These  predicted  amplitude  spectra  can  be  used  in  the  normaliza¬ 
tion  of  network  capability  simulations. 

•  Extrapolation  of  detection  capability  to  other  regions.  Since  we  have  an  estimate 
for  the  attenuation  for  paths  to  NORESS,  we  can  extrapolate  the  empirical  NORESS 
detection  capability  to  other  areas  where  the  attenuation  has  been  estimated. 

•  Regional  event  identification  Identification  of  regional  events  requires  that  we  dis¬ 
tinguish  relatively  small  differences  in  source  spectra.  This  requires  an  accurate, 
specific  knowledge  of  path  effects  between  source  and  array  locations. 


•  Yield  estimation.  If  regional  phases  recorded  at  NORESS  are  to  be  used  for  yield 
estimation,  their  attenuation  functions  must  be  accurately  determined. 

While  numerous  studies  of  Lg  attenuation  have  been  conducted,  comparatively  few  analyses 
for  regional  P  waves  exist.  However,  a  successful  parameterization  of  regional  P  wave  attenua¬ 
tion,  particularly  at  high  frequency,  has  important  implications  for  Q  in  the  lithosphere  and  for 
many  practical  issues  in  treaty  monitoring  seismology.  Analysis  of  Pn  is  difficult  because  it 
samples  a  much  smaller  fraction  of  the  focal  sphere  than  Lg,  so  its  amplitude  is  more  sensitive 
to  source  radiation  pattern,  focussing  and  defocussing,  and  scattering. 

Our  inversion  includes  Pn  and  Lg  spectra  from  190  regional  events  recorded  at  NORESS  at 
ranges  between  200  and  1300  km.  The  frequency  band  for  die  Lg  inversion  is  1  to  7  Hz. 
This  is  because  the  Lg  signal/noise  is  less  than  one  above  7  or  8  Hz  for  most  of  the  events 
beyond  800  km  range.  Also,  at  shorter  distances  the  Lg  spectra  are  probably  contaminated  by 
Sn  coda  at  high  frequency.  The  Pn  spectra  were  inverted  between  1  and  IS  Hz  (the  upper 
limit  of  adequate  signal/noise  for  the  events  in  the  data  base).  We  examined  die  dependence  of 
our  Q  estimation  on  azimuth,  assumptions  of  source  spectral  shape,  and  our  geometric  spread¬ 
ing  assumptions.  We  find  that  the  data  are  adequately  represented  by  a  simple  to2  source  and  a 
single,  frequency-dependent,  Q  model. 

1.2  Outline  of  the  report 

This  report  is  divided  into  8  sections,  including  this  introduction.  Section  2  is  a  brief  descrip¬ 
tion  of  regional  Pn  and  Lg  phases  and  reviews  previous  studies  of  their  attenuation.  Section  3 
describes  the  generalized  inversion  method,  including  all  the  simplifying  assumptions  used  in 
the  inversion.  The  dependence  of  our  results  on  these  assumptions  is  the  subject  of  Section  6.3. 
Section  4  and  Appendix  A  describe  the  data  base  and  the  signal  processing  procedures. 

The  primary  results  of  the  study  are  given  in  Section  5.  Lg  attenuation  is  discussed  in  Section 
5.1  and  Pn  attenuation  in  Section  5.2.  Because  the  Pn  and  Lg  inversions  are  done  separately, 
an  important  check  is  the  consistency  of  the  inverted  source  parameters  from  the  two  phases 
(Section  5.3).  In  Section  5.4  we  relate  the  moments  from  the  inversion  to  local  magnitude, 
which  is  important  for  predicting  range-dependent  spectra  for  a  given  magnitude.  In  Section 
5.5  we  discuss  the  relation  between  corner  frequency  from  the  inversion  and  source  size.  Most 
of  the  events  studied  have  M i  <  3.0  and  thus  have  high  comer  frequencies.  Therefore,  the  data 
do  not  clearly  resolve  source  corner  frequency,  but  we  note  that  our  results  are  consistent  with 


those  from  similar  studies.  Appendix  B  tabulates  the  inversion  parameters  for  each  event  stu¬ 
died  and  includes  plots  comparing  the  theoretical  spectra  and  observed  spectra. 

Error  analysis  is  the  subject  of  Section  6.  Both  formal  inversion  errors  and  estimated  parameter 
accuracy  are  discussed.  Section  6.3  is  devoted  to  analyzing  the  dependence  of  the  results  on 
our  parameterization.  In  particular,  the  source  spectrum  and  geometric  spreading  assumptioas 
are  varied,  and  the  resulting  change  in  model  parameters  is  monitored.  Section  6.4  discusses 
observed  systematic  variations  in  source  moment  versus  Lg  amplitude  for  different  mines. 
Section  6.5  illustrates  some  examples  of  spectra  that  are  not  well-modeled  by  the  inversion 
results.  Fortunately,  these  constitute  only  a  small  fraction  of  the  events  investigated. 

Section  7  demonstrates  one  of  the  most  important  applications  of  this  work,  the  simulation  of 
the  spectra  from  hypothetical  events.  In  particular,  predicted  range-dependent  Pn  and  Lg  spec¬ 
tra  for  fixed  magnitude  are  presented.  Finally,  Section  8  presents  our  main  conclusions  and 
summarizes  the  results  of  this  study. 


2.  REGIONAL  Pn  AND  Lg  PHASES 


Regional  seismograms  recorded  along  continental  paths  are  characterized  by  the  appearance  of 
one  or  more  of  the  phases  Pn,  Pg,  Sn,  and  Lg.  Pn  and  Sn  predominantly  sample  the  upper¬ 
most  mantle  while  Pg  and  Lg  are  crustal  phases.  In  this  study  we  focus  our  attention  on  Pn 
and  Lg.  Lg  is  often  the  largest  amplitude  signal  for  stable  continental  paths.  It  has  been 
modeled  successfully  as  a  sum  of  higher  mode  surface  waves  with  constant  group  velocity  near 
3.5  km/s  and  energy  density  confined  to  the  crust  (Knopoff,  et  al.,  1973;  Panza  and  Calcagnile, 
1975;  Bache,  et  al.,  1981).  Regional  Pn  phases  are  the  first  arrivals  beyond  about  200  km  and 
propagate  in  the  upper  mantle  with  group  velocities  typically  between  7.5  and  8.0  km/s.  How¬ 
ever,  at  distances  less  than  300  km,  Pg  arrives  within  5s  of  Pn.  The  Pn  amplitude  is  gen¬ 
erally  less  than  that  of  Lg  (by  a  factor  of  3  or  more  in  the  distance  range  100-500  km;  Myk- 
keltveit  and  Ringdal,  1979),  due  partly  to  geometric  spreading  differences,  and  it  has  a  higher 
dominant  frequency  than  Lg. 

Numerous  studies  of  Lg  attenuation  have  been  conducted  (e.g.,  Cheng  and  Mitchell,  1981; 
Herrmann  and  Kijko,  1983;  Singh  and  Herrmann,  1983;  Peseckis  and  Pomeroy,  1984;  Cam- 
pillo,  et  al.,  1985;  Hasegawa,  1985;  Chun,  et  al.,  1987;  Gupta  and  McLaughlin,  1987;  Shin  and 
Herrmann,  1987).  Table  2.1  and  Figure  2.1  summarize  reported  Lg  attenuation  estimates  for 
various  regions.  Entries  1-4  of  Table  2.1  are  attenuation  estimates  for  paths  in  the  western 
United  States.  With  the  exception  of  Chavez  and  Priestley  (1986),  these  studies  are  band  lim¬ 
ited  to  relatively  low  frequencies.  Entries  5-10  are  estimates  for  paths  across  eastern  North 
America.  An  obvious  conclusion  is  that  the  tectonically  active  western  United  States  is  charac¬ 
terized  by  a  lower  Q  and  stronger  frequency  dependence  than  the  stable  eastern  North  Ameri¬ 
can  shield. 

Very  few  regional  P  wave  attenuation  studies  have  been  conducted.  This  is  probably  because 
Pn  samples  a  much  smaller  fraction  of  the  focal  sphere  than  Lg,  making  its  amplitude  more 
sensitive  to  source  radiation  pattern,  focussing  and  defocussing,  and  scattering.  Several  recent 
studies  of  both  P  and  S  wave  attenuation  have  uncovered  the  peculiar  observation  that  Qp/Q0 
reaches  values  greater  than  or  equal  to  one  at  high  frequency  (Clements,  1982;  Taylor,  et  al., 
1986;  Butler,  et  al.,  1987;  Chavez  and  Priestley,  1987;  Sereno  and  Orcutt,  1987).  In  a  study  of 
regional  crustal  phases,  Chavez  and  Priestley  (1987)  estimated  attenuation  between  I  and  10 
Hz  for  paths  across  the  Great  Basin  and  found  that  Pg  Q  was  less  than  Lg  Q.  We  have  exam¬ 
ined  both  Pn  and  Lg  attenuation  in  this  study,  but  note  that  the  two  phases  sample  different 
depths.  Both  the  scattering  and  intrinsic  absorption  contributions  to  the  apparent  attenuation 


Table  2.1.  Reported  Qi^f)  estimates. 


Region 

Frequency  (Hz) 

Ql,U) 

Co  n 

Reference 

(1) 

Western  United  States 

0.5-3.S 

140-200 

0.4-0.6 

Singh  and  Herrmann  (1983) 

(2) 

Western  United  States 
(NTS  Explosions) 

1. 0-2.0 

200-300 

0.2-0.4 

Peseckis  and  Pomeroy  (1984) 

(3) 

Western  United  States 
(NTS  Explosions) 

1. 0-2.0 

139 

0.6 

Nuttli  (1986) 

(4) 

Western  United  States 
(Great  Basin) 

0.25-12.5 

206 

0.68 

Chavez  and  Priestley  (1986) 

(5) 

Eastern  United  States 

0.5-3.5 

1000 

0.3-0.4 

Singh  and  Herrmann  (1983) 

(6) 

Eastern  North  America 

1.0-15.0 

900 

0.2 

Hasegawa  (1985) 

(7) 

Eastern  United  States 

— 

1000 

0.35 

Goncz  and  Dean  (1986) 

(8) 

Eastern  Canada 

0.6-10.0 

1100 

0.19 

Chun,  er  al.  (1987) 

(9) 

Eastern  United  States 

0.5-7.0 

800 

0.32 

Gupta  and  McLaughlin  (1987) 

(10) 

Eastern  North  America 

0.5-15.0 

500-550 

0.65 

Shin  and  Herrmann  (1987) 

(11) 

Central  France 

0.5-10.0 

290 

0.52 

Campillo,  et  al.  (1985) 

(12) 

Southern  Africa 

— 

600 

0.4 

Mitchell,  et  al.  (1987) 

are  expected  to  be  different  for  the  two  phases.  Therefore,  comparison  of  the  Q  obtained  does 
not  provide  much  information  about  lithosphere  rheology.  Numerical  modeling  of  Lg  suggests 
that  its  observed  attenuation  is  an  effective  measure  of  the  average  absorption  of  shear  waves 
in  the  crust  (Campillo,  et  al.,  1985).  Factors  other  than  intrinsic  absorption  play  an  important 
role  in  the  attenuation  of  Pn.  Therefore,  it  is  inappropriate  to  associate  observed  Pn  attenua¬ 
tion  with  intrinsic  absorption  at  a  particular  depth  in  the  upper  mantle.  Our  emphasis  is  not, 
however,  on  the  physical  mechanisms  of  attenuation,  but  rather  on  the  development  of  empiri¬ 
cal  functions  to  be  used  to  predict  frequency-dependent  amplitudes. 

Regional  Lg  phases  have  also  been  used  to  estimate  source  parameters  including  comer  fre¬ 
quency  and  seismic  moment  (Street,  et  al.,  1975;  Dwyer,  et  al.,  1983;  Hasegawa,  1983;  Shin 
and  Herrmann,  1987).  These  results  are,  in  general,  consistent  with  near-field  source  studies, 
with  the  possible  exception  that  corner  frequencies  obtained  from  Lg  tend  to  be  lower  than 
those  estimated  at  short  range  (Mueller  and  Cranswick,  1985). 

Observed  Lg  to  P  amplitude  ratios  are  complicated  functions  of  source  depth,  near-surface 
velocity,  and  lateral  structural  variations.  Nevertheless,  moderate  success  has  been  attained  in 
using  LgIPn  and  Lg/Pg  ratios  as  regional  earthquake-explosion  discriminants  (for  review,  see 
Pomeroy,  et  al.,  1982).  For  example,  in  a  comparison  of  LgIPn  ratios  for  an  NTS  explosion 
to  a  co-located  earthquake  at  a  range  of  450  km,  Willis  (1963)  found  that  the  earthquake 
source  resulted  in  an  LgIPn  ratio  5  times  that  of  the  explosion.  Pomeroy  (1977)  found  a  simi¬ 
lar  result  comparing  12  earthquakes  to  the  SALMON  nuclear  explosion  detonated  in  Missis¬ 
sippi.  More  recently,  Murphy  and  Bennett  (1982)  compared  NTS  explosions  to  nearby  earth¬ 
quakes  and  found  that  observed,  Lg  amplitudes  were  typically  greater  for  earthquakes  than  for 
explosions  with  comparable  P  wave  amplitudes.  However,  they  found  that  the  simple  time 
domain  amplitude  ratios  did  not  provide  consistent  separation  between  the  two  populations. 
Explosions  generate  Lg  energy  primarily  through  P-SV  mode  conversions  and  scattering.  In 
addition,  earthquakes  directly  produce  shear  wave  energy  that  can  contribute  to  Lg.  Therefore, 
for  a  given  focal  depth  earthquakes  are  expected  to  be  more  efficient  in  Lg  excitation  than 
explosions.  Bennett,  et  al.  (1987)  compared  the  relative  excitation  of  synthetic  Pg  and  Lg 
phases  for  earthquakes  and  explosions  as  a  function  of  source  depth.  They  found  that  the 
Lg/Pg  ratio  was  only  slightly  higher  for  near-surface  earthquakes  than  for  near-surface  explo¬ 
sions,  but  was  a  factor  of  3  higher  at  a  depth  of  1  km. 


3.  GENERALIZED  INVERSION 


The  key  to  estimating  regional  wave  attenuation  is  separating  source  and  path  effects.  For  one 
event  recorded  at  one  station,  there  is  a  direct  trade-off  between  source  and  path  contributioas 
to  the  observed  signal  spectrum.  We  have  many  events  recorded  at  one  station  and  study  the 
attenuation  as  a  function  of  distance  and  frequency.  Our  analysis  assumes  that  all  observed 
spectra  can  be  fit  by  a  single  frequency -dependent  Q  model  and  that  all  source  spectra  scale 
simply  with  event  size.  In  this  way  we  use  the  range-independence  of  the  source  function  to 
separate  the  different  contributions  to  the  seismic  spectrum. 

Our  generalized  inversion  procedure  simultaneously  estimates  seismic  moment  and  Q(f). 
Adopting  standard  methods  for  solving  non-linear  inverse  problems,  we  linearize  the  system  of 
equations  governing  the  relationship  between  the  data  and  model  parameters.  We  assume  a 
starting  model,  compute  theoretical  data,  subtract  it  from  the  observed  data,  and  solve  itera¬ 
tively  for  the  model  perturbations  that  minimize  the  data  residual  in  the  least-squares  sense.  In 
practice,  we  have  found  it  necessary  to  include  damping  to  stablize  the  solution.  That  is,  we 
minimize  a  weighted  sum  of  the  data  residuals  and  the  model  perturbation  norm. 


3.1  Description  of  the  method 

We  parameterize  the  instrument-corrected  amplitude  spectrum  of  a  seismic  signal  as 


A(f,r)  =  S(J)  G(r,r0)  exp 


(3.1.1) 


where  A(f,r)  is  the  observed  displacement  spectrum  at  range  r  and  frequency  /;  S(f)  is  the 
source  spectrum;  G(r,ro)  is  geometric  spreading,  and  the  last  term  is  the  effective  attenuation 
for  travel  time,  t.  Note  that  the  effective  attenuation  includes  contributions  from  both  anelasti- 
city  and  scattering. 


We  assume  the  geometric  spreading  function,  so  the  computed  attenuation  functions  and 
seismic  moments  are  relative  to  the  assumed  spreading  rate.  Following  Herrmann  and  Kijko 
(1983),  we  express  the  spreading  function  as 

G(r,r0)  =  (i/r)>  for  r  <  r0 

=  r~Q  (ro/r)m,  for  r  >  r0  (3  ,  2) 


8 


where  r0  is  a  transition  distance  from  spherical  spreading  to  spreading  rate  m.  By  comparing 
the  long  period  amplitude  spectrum  of  Lg  to  moments  calculated  from  long  period  surface 
waves,  Street,  et  al.  (1975)  empirically  determined  r0  =  100  km,  or  roughly  twice  the  crustal 
thickness.  Measuring  the  decay  rate  of  synthetic  Lg  phases  computed  for  an  elastic  medium, 
Herrmann  and  Kijko  (1983)  verified  that  Lg  frequency  domain  spreading  was  accurately 
described  as  cylindrical  (m=l/2)  and  substantiated  the  empirical  result  of  Street,  et  al.  (1975) 
for  r0.  Less  work  has  been  done  on  the  spreading  rate  of  Pn.  Because  its  energy  density  is 
more  localized  about  a  single  ray  path,  Pn  geometric  spreading  is  more  sensitive  to  velocity 
gradients  in  the  upper  mantle.  Numerical  studies  of  Pn  indicate,  for  typical  upper  mantle 
structures,  that  its  spreading  rate  is  lies  between  r-1  and  r~~  (Langston,  1982;  Wallace,  per¬ 
sonal  communication).  This  approximately  corresponds  to  the  range  between  simple  turning 
ray  and  canonical  headwave  Pn  ray  path  interpretations.  In  this  report,  we  investigate  spread¬ 
ing  rates  between  these  bounding  values  as  approximations  to  the  true  spreading  rate  of  Pn. 


The  inversion  also  requires  a  source  spectral  shape  parameterization.  We  assume  an  or  spec¬ 
trum  with  corner  frequency  inversely  proportional  to  the  cube-root  of  the  long  period  source 
level.  That  is,  we  assume  the  seismic  source  function  can  be  expressed  as 


S(f)  =  Sq  s(f ) 


(3.1.3) 


where  50  is  the  long  period  source  level  and  s(J)  is  a  known  function  which  depends  on  50  and 
describes  the  source  spectral  shape.  We  have  investigated  two  functional  forms  for  the  to2 
source;  a  simplified  Mueller-Murphy  (1971)  explosion  source  model  and  a  Brune  (1970,  1971) 
earthquake  source  model.  The  explosion  source  has  the  form 


*(/)  = 


fc 


(fc4  +  (1-2(3) /c  2/2  +  P2/4)1' 


(3.1.4) 


where  (}  controls  the  amount  of  overshoot  and  fc  is  the  corner  frequency.  The  earthquake 
source  model  has  frequency  dependence 


s{f)  = 


fc2 

(f  +  fc)2 


(3.1.5) 


In  both  cases  we  assume  that  corner  frequency  scales  with  the  cube-root  of  the  long  period 
levei.  That  is, 


fc  =  cS0 


-1/3 


(3.1.6) 
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and  c  is  a  parameter  of  the  inversion.  In  Figure  3.1  we  plot  a  family  of  curves  for  each  source 
model. 

The  source  parameters  estimated  by  the  inversion  are  Sq  for  each  event  and  a  single  value  of  c 
relating  corner  frequency  to  long  period  source  level.  For  near-surface  explosions,  the  relation¬ 
ship  between  Sq  and  seismic  moment  ( Mq )  for  Pn  is  (Stevens  and  Day,  1985) 

Moxp 

S>QXp(Pn)  = - -  —  (3.1.7) 

47TpX 

All  of  the  explosions  in  our  data  set  are  near-surface  mining  explosions,  characterized  by  the 
near-surface  density  (p5)  and  compressional  velocity  (a*).  For  earthquakes  the  expression  is 
complicated  by  the  depth  of  the  event  and  the  radiation  pattern.  From  Stevens  and  Day  (1985), 
the  expression  for  Sq  is 


SfftPn)  = 


4n(PcPsacas)l/2 


(3.1.8) 


where  pt.  and  ac  are  crustal  density  and  compressional  velocity  at  the  source  depth  and  Rq0  is 
the  P  wave  source  radiation  pattern.  Since  our  observations  are  from  a  single  station,  and  the 
focal  mechanisms  for  the  small  earthquakes  in  our  data  set  are  unknown,  we  do  not  know  R^0 
and  are  therefore  unable  to  estimate  seismic  moment  from  S^iPn).  However,  using  the  empiri¬ 
cal  result  of  Street,  et  al.  (1975),  we  can  estimate  the  earthquake  moments  from  Lg  spectra 
using 


W-£)  = - -r,  (3.1.9) 

4*PcPc 

where  pc  and  are  the  average  crustal  density  and  shear  wave  velocity,  respectively. 

The  amount  of  Lg  energy  excited  by  an  explosion  is  depth  dependent  and  complicated  by 
near-source  wave  conversions  Also,  a  cviindrically  symmetric  surface  explosion  generates 
relatively  more  shear  wave  energy  than  one  that  is  buried.  Thus,  a  simple  expression  relating 
the  seismic  moment  of  an  explosion  to  the  long  period  Lg  spectrum  does  not  exist.  In  general, 
earthquakes  generate  more  shear  wave  energy  than  explosions,  so  we  express  the  long  period 
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Lg  source  level  for  an  explosion  as  an  unknown  fraction  of  the  long  period  level  for  an  earth¬ 
quake  of  equal  moment  and  equal  depth.  That  is, 

StfHLg)  =  k  Stf(Lg)  (3.1.10) 

where  K  is  an  unknown  constant,  presumably  less  than  one.  We  will  estimate  this  constant  by 
using  the  explosion  moments  obtained  from  Pn  in  the  above  expression  for  SQKp(Lg). 

The  inversion  requires  a  parameterization  of  apparent  attenuation.  No  attempt  is  made  to 
separate  the  effects  of  scattering  from  intrinsic  absorption.  We  simply  model  the  spectral  decay 
in  terms  of  a  power-law  frequency  dependence  of  Q, 

Q(f)  =  QofA  (3.1.11) 

with  ()0  and  r|  parameters  of  the  inversion.  Using  (3.1.1)  and  (3.1.1 1)  we  compute  the  function 
d(f,r)  from 


d{f,r)  =  Log  A(f,r)  -  Log  G(r,r0 )  =  Log  S(J)  -  O.Ola0/(1  T|)  t  (3.1.12) 


where 


100  K  Log  e  , 

«o  =  - .  *  (3.1.13) 

C/o 

The  d(f,r)  is  the  observed  data  corrected  for  the  assumed  geometric  spreading  function  and  is 
the  input  data  for  the  inversion.  The  factor  of  100  in  (Xq  was  introduced  to  avoid  matrix  ill- 
conditioning.  A  d^(f,r)  can  be  computed  from  an  assumed  starting  model.  This  is  subtracted 
from  the  observed  data  d(f,r),  giving  a  data  residual  which  can  be  expressed  in  matrix  form  as 

Ad  =  A  Am  (3.1.14) 

In  (3.1.14),  Ad  is  the  data  residual  vector,  A  is  the  matrix  of  partial  derivatives  of  the  data 
with  respect  to  the  model  parameters,  and  Am  is  the  model  perturbation  vector.  Explicitly 


Ad1  -  [Ad{fx,rx),  AJ(/2,r, ),...,  Adifnf,rx ),  Ad(j\,r2),  Ad(f2,r2) . Ad(fnj,rnt.v)\ 


(3.1.15) 


where  nf  is  the  number  of  frequencies  and  nev  is  the  number  of  events.  The  total  number  of 
data  is  nd  =  nf  x  nev  and  the  number  of  parameters  is  up  =  nev  +  3.  The  {np  x  1)  model 
perturbation  vector  is 

Am7  =  [A (LogSl0),  A(LogSl),...,  A {LogStfv),  Ac,  Aa0,  An)  (3. 1.16) 

where  Sq  is  the  long  period  source  level  of  the  i,h  event.  The  first  nev  columns  of  the 
(nd  x  np)  matrix  A  are  the  partial  derivatives  of  the  data  with  respect  to  LogSy  tor  each  event. 
The  last  3  columns  of  A  are  the  partials  with  respect  to  the  parameters  c,  (Xo,  and  ry  That  is, 

Bd^  ddi  ddi  Bd^  Bd\  Bd^ 

BLogSl0  BLogSl  '  '  dLogStfv  dc  Ba0  9n 

A=  .  .  .  (3.1.17) 

Bd^  BdrlJ 
BLogSfQV  3c  Ba0  B  n 

Each  of  the  first  nev  columns  of  A  has  only  nf  non-zero  elements.  That  is.  spectra  from  the  i,h 
event  have  no  dependence  on  the  long  period  level  of  the  jlh  source,  unless  i  =  j  That  part  of 
the  matrix  A  is  block  diagonal. 


BLogS'0  BLogSl 


The  partials  with  respect  to  path  parameters  are 


=  -0.01  /(1~T1)  r. 


BaQ 


Bd 


~  =-0.01  aoMln/)/(1-^. 

dri 


( 3.1.18) 

(3 .1.19) 


The  absolute  value  of  Bd/B(Xo  typically  ranges  between  0.6  and  8  for  Lg  and  between  0.3  and 
5  for  Pn.  The  partial,  Bd/Bt\,  varies  between  0  and  7  for  Lg  and  0  and  6  for  Pn 


For  the  Mueller-Murphy  (1971)  source  model,  the  partials  with  respect  to  source  terms  are 


Bd 


BLogS0  3 


2(52/45^  +  (1-2P)c*:/‘5B 


2f  2 o-l/3 


c4  +  (l-2B)c2/  3  +  P2/45g/3  ' 


1 3  1.20) 


dd  _  2  Log  e 
3c  c 


2c3  +  (l-2ft)c/2S5/3 

+  (l-2P)c2/25^  +  p2/4.Vf)> 


0.1.2!) 


Similarly,  for  the  Brune  tl970,  1971)  source  model 

SJ  =  t 

DLofjio  3(Sj'3/+  «)  ’ 

■  ■ 
dd  _  2  Lo#  e _ c 

dc  c  Sq*  +  c 

The  partials  with  respect  to  LogSQ  are  equal  to  one  for/ C  fc  and  1/3  for / >•  /..  The  partials 

with  respect  to  c  are  zero  at  zero  frequency  and  increase  to  a  maximum  value  of  for 

c 

f  :*»  fc.  Examples  of  partials  with  respect  to  the  two  source  parameters  are  displayed  in  Figure 
3.2. 


(3.1.22) 

(3.1.23) 


The  system  of  equations  defined  by  (3.1.14)  is  over-determined,  and  an  exact  solution  is  not 
available.  The  solution  that  minimizes  the  data  residual  in  the  least-squares  sense  is  given  by 
the  normal  equations 


Am  =  (ATAf]ATAd  (3.1.24) 

The  damped  least-squares  solution,  for  damping  factor  X,  is 

Am  =  (ArA  +  >JflATAd  (3.1.25) 

The  actual  matrix  inversion  was  done  using  the  boardering  method  (Press,  et  al.  1986).  Singu¬ 
lar  value  decomposition  was  also  used  to  solve  the  generalized  inversion,  but  was  found  to 
give  the  same  result  as  boardering,  which  can  be  done  much  faster. 

3.2  Assumptions 

In  this  section  we  discuss  the  simplifying  assumptions  used  in  the  inversion  The  dependence 
of  the  results  on  these  assumptions  is  investigated  later  in  Section  6  3. 
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3.2.1  Source  spectra 

We  assume  that  the  source  spectrum  is  uniquely  defined  by  its  long  period  level.  This  is  a 
common  assumption  which  has  a  controlling  influence  in  some  studies  of  Q.  However,  it  is  not 
crucial  for  our  analysis  because  most  of  the  events  studied  have  Mi  <  3.0,  so  their  corner  fre¬ 
quencies  are  near  or  beyond  the  upper  limit  of  bandwidth  inverted.  For  example,  the  bandwidth 
for  Lg  is  1-7  Hz,  and  we  found  that  inversion  with  a  flat  source  spectrum  gives  nearly  the 
same  Q  as  one  with  the  more  realistic  source  functions. 

For  the  characteristics  of  the  source  beyond  the  corner  frequency,  we  follow  the  practice  of  the 

•y 

vast  majority  of  the  seismological  literature  and  assume  the  source  spectrum  follows  an  or 
decay.  The  theoretical  basis  for  this  assumption  is  not  very  strong,  but  careful  empirical  stu¬ 
dies  generally  support  its  validity.  For  example,  Chael  (1987)  studied  12  aftershocks 
(3.3  <  WbLg  ^  5.8)  of  the  1982  Miramichi  earthquake  and  found  that  the  data  strongly  favor 
or  over  (d3 * 5  source  models.  He  used  spectral  ratios  of  large  to  small  co-located  events  to  elim¬ 
inate  the  shared  path  contribution  and  to  measure  the  spectral  decay  between  the  two  corner 
frequencies.  Another  issue  is  the  scaling  of  corner  frequency  with  event  size,  which  is  a  sub¬ 
ject  of  numerous  studies  in  the  literature.  The  most  straightforward  assumption  is  cube-root 
scaling,  though  it  is  recognized  that  the  scaling  can  also  be  by  a  smaller  exponent  due  to 
changes  in  source  depth  (explosions)  or  stress  drop  (earthquakes)  that  correlate  with  event  size. 

Specifically,  we  have  incorporated  the  Brune  (1970,  1971)  earthquake  source  model  (3.1.5)  and 
the  Mueller-Murphy  (1971)  explosion  model  (3.1.4)  into  the  inversion.  Most  of  the  events  in 
our  data  base  are  small,  near-surface  mining  explosions  that  appear  to  have  corner  frequencies 
higher  than  7  Hz.  In  Section  5,  results  for  the  Mueller-Murphy  model  are  presented  and  in 
Section  6.3  we  study  the  effect  of  altering  our  source  assumptions. 

3  2  2  Geometric  spreading 

The  geometric  spreading  of  regional  Lg  phases  is  well-constrained,  both  empirically  and 
theoretically  The  onset  consists  of  higher  mode  surface  waves  which  are  accurately  described 
at  long  ranges  by  cylindrical  spreading.  The  choice  of  r0  in  (3.1.2),  the  transition  distance 
between  spherical  and  cylindrical  spreading,  was  estimated  empirically  by  Street,  et  al.  (1975) 
It  was  later  substantiated  theoretically  by  Herrmann  and  Kijko  (1983)  as  approximately  twice 
the  crustal  thickness  We  adopt  their  value  of  r{)  =  KX)  km  lor  the  results  presented  in  Section 

5  Because  all  of  the  events  used  this  study  were  at  ranges  greater  than  r0,  its  value  has  no 
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effect  on  the  Q  estimate,  but  it  trades  off  directly  with  the  seismic  moment  estimate  for  each 
event. 

The  spreading  rate  of  regional  Pn  phases  is  not  well-constrained.  Pn  samples  a  small  fraction 
of  the  focal  sphere,  and  its  geometric  spreading  rate  is  dependent  upon  velocity  gradients  in  the 
upper  mantle.  An  important  constraint  on  Pn  spreading  is  that  the  ratio  between  Lg  and  Pn 
long  period  source  levels  be  range-independent  This  criterion  supports  a  choice  of  r-1  ?  for  the 
Pn  spreading  rate,  and  this  is  used  for  the  results  of  Section  5.  Results  for  other  spreading 
rates  between  r_l  and  r~ ^  are  investigated  in  Section  6.3. 

3.2.3  Q(f)  parameterization 

We  choose  a  simple  parameterization  of  the  range-dependent  decay  of  the  seismic  spectrum  in 
terms  of  a  power-law  frequency  dependence  of  Q.  We  do  not  attempt  to  distinguish  intrinsic 
absorption  from  scattering.  That  is,  we  explicitly  acknowledge  that  our  Q(f)  is  an  empirical 
parameterization  of  the  data  and  do  not  address  the  interpretation  in  terms  of  rheology.  The 
justification  of  this  parameterization  is  its  compatibility  with  other  studies  (for  comparison)  and 
the  fact  that  it  does  allow  us  to  obtain  satisfactory  agreement  between  observed  and  theoretical 
spectra. 

3  2.4  Radiation  pattern 

In  assuming  that  corner  frequency  is  inversely  proportional  to  the  cube-root  of  the  long  period 
source  level,  we  have  assumed  isotropic  radiation.  While  this  is  reasonable  for  explosions,  it 
could  result  in  significant  error  for  earthquakes.  A  large  event  near  a  P- wave  nodal  plane  will 
have  a  small  SQ{Pn)  and  a  low  corner  frequency,  which  cannot  be  modeled  with  our  parame¬ 
terization.  Focal  mechanisms  are  not  available  for  the  small  events  in  our  data  set.  We  have 
run  the  inversion  on  a  subset  of  known  explosions,  and  on  the  full  data  set,  with  most  of  the 
added  events  probably  being  earthquakes.  We  note  that  the  variance  of  the  inversion  increases 
for  the  full  data  set,  presumably  due  to  uninodeled  radiation  pattern  effects.  However,  the 
results  do  not  change  significantly,  which  indicates  that  the  radiation  pattern  is  not  biasing  the 
estimated  parameters.  Because  Lg  samples  a  larger  fraction  of  the  focal  sphere,  the  radiation 
pattern  is  of  less  concern  for  Lg  than  for  Pn. 

3  2.5  Interference  phenomena 

The  constructive  and  destructive  interference  of  multiple  arrivals  can  modulate  the  observed 
seismic  spectrum.  Examples  include  spectral  scalloping  resulting  from  interfering  depth  phases. 
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local  site  resonances,  and  multiple  sources.  Time  lags  greater  than  approximately  0.5s  are  of 
little  concern,  since  frequency  smoothing  suppresses  that  effect.  Multiples  with  lags  less  than 
150  ms  can,  however,  degrade  the  accuracy  of  the  Q  estimate.  As  an  example,  Baumgardt  and 
Ziegler  (1987)  attributed  a  broad  spectral  peak  centered  near  7  Hz  to  ripple-firing  of  mine 
blasts  with  delays  of  the  order  of  150  ms.  In  this  case,  destructive  interference  could  increase 
the  observed  spectral  decay  rate  between  approximately  7  and  14  Hz.  However,  for  this  to 
have  a  significant  impact  on  our  results,  ripple-firing  with  very  coexistent  time  lags  at  many 
different  mine  locations  would  have  to  occur.  For  time  lags  greater  than  150  ms,  the  inversion 
finds  a  smooth  curve  that  fits  through  the  modulated  spectra.  Many  examples  are  pictured  in 
Appendix  B  . 


4.  DATA  BASE 


The  data  used  in  this  study  consist  of  stable,  array-averaged  spectra  for  190  regional  events 
recorded  by  the  small  aperture  NORESS  seismic  array  in  Norway.  The  NORESS  array 
configuration  and  sampling  rate  were  designed  to  enhance  the  detection  of  signals  from  small 
regional  events  (Mykkeltveit,  et  al.t  1983).  The  array  includes  25  short  period  instruments  in 
concentric  rings  with  maximum  diameter  of  3  km.  The  data  are  digitally  recorded  at  40 
samples/s.  Figure  4.1  shows  the  NORESS  array  configuration  and  the  short  period  instrument 
response  to  die  Nyquist  frequency. 

4.1  Signal  processing 

The  calculation  of  seismic  spectra  is  incorporated  into  an  automated  seismic  array  processing 
program  ( SAIAP )  developed  at  SAIC  as  an  extension  of  the  RONAPP  program  used  at  NOR- 
SAR  (Mykkeltveit  and  Bungum,  1984).  The  program  computes  spectra  for  each  automatically 
detected  signal.  The  spectral  estimation  technique  is  that  proposed  by  Bache,  et  al.  (1985).  A 
10%  cosine-squared  taper  is  applied  to  a  5s  window  starting  0.3s  before  the  onset  time  of  the 
arrival  on  the  vertical  component.  The  time  series  is  padded  with  zeros  to  1024  samples  and 
fast  Fourier  transformed.  The  same  procedure  is  applied  to  a  noise  sample  taken  prior  to  the 
first  P  detection.  The  squared  noise  amplitude  spectrum  (power)  is  subtracted  from  the  squared 
signal  spectrum  (energy  density).  The  resulting  noise-corrected  signal  spectra  are  averaged 
across  the  array  and  corrected  for  the  instrument  response.  Bache,  et  al.  (1985)  demonstrated 
that  if  the  noise  is  random,  stationary,  and  uncorrelated  with  the  signal,  the  signal  spectrum 
estimate  obtained  with  this  method  converges  to  the  true  signal  spectrum  as  the  number  of  ele¬ 
ments  increases.  Array-averaging  has  the  desirable  effect  of  suppressing  uncorrelated  near¬ 
receiver  local  site  effects. 

Since  the  seismic  spectra  are  computed  with  fixed  time  window  lengths,  they  encompass 
different  group  velocity  windows  at  different  ranges.  We  note  that  the  spectra  computed  this 
way  may  be  different  from  those  computed  with  fixed  group  velocity  windows,  particularly  for 
Lg,  where  short  range  spectra  contain  contributions  from  more  modes  than  long  range  spectra 
which  are  confined  to  a  smaller  group  velocity  window.  For  example,  Chun,  et  al.  (1987) 
estimated  vertical  component  Lg  attenuation  using  both  fixed  window  lengths  and  fixed  group 
velocity  windows  on  the  same  data  set.  Using  ECTN  (Eastern  Canada  Telemetered  Network) 
Lg  data,  they  found  Q(J)  =  800/ 026  for  a  fixed  time  window  of  17.07s.  For  a  fixed  group 
velocity  window  of  3.61  to  2.6  km/s,  they  found  Q(f)  =  1 100/°  19  Their  data  spanned  dis¬ 
tances  of  90  to  867  km,  so  the  time  windows  varied  from  9.7s  to  93s.  It  is  not  clear  whether 
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the  fixed  time  window  gives  higher  attenuation  because  the  lower  order,  low  group  velocity 
modes  are  not  included  at  large  distance  (the  explanation  preferred  by  Chun,  et  al.,  1987)  or 
because  increased  scattered  energy  and/or  greater  noise  contamination  is  included  in  the  long 
time  segments.  The  situation  is  further  complicated  by  the  increased  possibility  that  .spurious 
arrivals  are  included  in  the  long  window.  In  any  case,  we  have  chosen  fixed  time  windows 
because  they  are  most  convenient  for  automatic  processing,  and  our  results  correctly  represent 
the  spectral  character  of  Lg  windowed  in  this  way. 

As  a  final  step  in  the  signal  processing,  the  log  amplitude  spectra  are  smoothed  over  a  2  Hz 
frequency  band.  As  an  example,  Figure  4.2  compares  the  smoothed  Pn  spectrum  to  the  array- 
averaged  spectrum  obtained  by  the  automatic  processing.  The  dashed  curve  is  the  average 
noise  spectrum  estimated  from  samples  taken  prior  to  Pn  for  many  events  (Henson  and  Bache, 
1986). 

4.2  Data 

Table  4.1  is  a  list  of  all  events  and  phases  (Pn  and/or  Lg)  used  in  the  inversion.  In  some  cases 
only  one  phase  was  included  because  the  other  was  not  detected  or  because  it  had  a  low 
signal/noise  ratio  over  the  frequency  band  used  in  the  inversion.  Reported  mining  explosions 
are  identified  by  an  "EX"  under  the  column  heading,  TYPE.  The  label  following  "EX" 
identifies  the  mine.  Three  are  in  southwest  Norway;  BLA  (Blasjo),  TIT  (Titania),  and  NYG 
(Nygardstaugen).  Other  mine  codes  are  those  used  in  the  bulletin  published  by  the  University 
of  Helsinki  based  on  the  Finish  Seismic  Array.  The  mine  locations  are  listed  in  Table  4.2, 
along  with  distance  and  azimuth  from  NORESS,  and  are  displayed  in  Figure  4.3.  An  "EQ" 
designator  identifies  presumed  earthquakes,  although  some  may  be  unreported  explosions. 
Events  that  are  not  reported  explosions,  but  have  locations  within  50  km  of  known  mines  are 
considered  of  unknown  source  type.  The  only  exception  is  event  10  of  Table  4.1,  which  was 
labeled  by  a  NORSAR  analyst  as  a  probable  dam  explosion. 

The  location  and  origin  times  are  from  a  local  bulletin  published  by  the  University  of  Bergen 
or  the  University  of  Helsinki,  when  available,  or  from  the  Preliminary  Determination  of  Epi¬ 
centers  (PDE)  bulletin.  Events  for  which  an  independent  network  solution  is  not  available  are 
assigned  either  RONAPP  (single  asterisk)  or  SAIAP  (double  asterisk)  locations.  The  RONAPP 
locations  are  published  in  the  NORESS  bulletin,  along  with  the  detection  times  of  the  P  and 
Lg  phases  used  in  the  location  solutions.  In  some  cases  RONAPP  did  not  associate  phases 
correctly,  and  we  were  able  to  use  the  SAIAP  solution  with  the  appropriate  phase  association. 
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Figure  4.2.  Amy-averaged  displacement  spectrum  of  a  ML  -  3.2  Estonia  mine  blast  at  a 
range  of  930  km.  Superimposed  is  the  same  spectrum  smoothed  over  a  2  Hz 
bandwidth.  The  dashed  curve  represents  average  NORESS  noise  estimated  from 
samples  taken  prior  to  Pn  for  many  events  (Henson  and  Bache,  1986). 
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Table  4.1.  Events  used  in  the  generalized  inversion. 

£  Event 

Date  Time  Location  Type  Phases 

10-25-85 

10-27-85 

10-27-85 

10-27-85 

10-29-85 

10-31-85 

10- 31-85 

11-  6-85 
11-  9-85 
11-  9-85 
11-12-85 
11-13-85 
11-13-85 
11-14-85 
11-15-85 
11-21-85 
11-21-85 
11-23-85 
11-25-85 
11-27-85 
11-27-85 
11-28-85 

11- 30-85 

12-  1-85 
12-  5-85 
12-  7-85 
12-  7-85 
12-  7-85 
12-10-85 
12-10-85 
12-10-85 
12-11-85 
12-11-85 
12-13-85 
12-14-85 
12-17-85 
12-23-85 
12-24-85 
12-24-85 
12-25-85 
12-25-85 
12-25-85 
12-27-85 
12-27-85 
12-27-85 


2:04 
4:36 
4:41 
4:52 
0:23 
2:56 
4:11 
4:51 
4:43 
8:21 
2:22 
2:08 
4:11 
2:52 
13:54 
11:50 
13:17 
13:06 
13:06 
4:54 
12:18 
9:30 
19:05 
7:21 
12:25 
13:18 
14:16 
14:39 
12:06 
12:18 
13:43 
12:14 
12:51 
12:09 
14:35 
13:08 
4:27 
0:04 
13:13 
12:04 
13:19 
14:18 
11:06 
12:16 
12:42 


59.30N  28.10E 
61.30N  4.30E 
66.40N  1 1.60E* 
66.00N  14.10E** 
59.31N  6.95E 
62.78N  18.03E 
60.70N  29.00E 
59.31N  6.95E 
58.34N  6.43E 
62.60N  6.70E** 
59.50N  25.00E 
59.30N  28.10E 
58.34N  6.43E 
60.70N  28.70E 
61. ION  29.90E 
59.30N  27.20E 
59.29N  7.04E 
59.50N  25.00E 
59.40N  28.50E 
59.73N  5.71E 
61.40N  31.60E 
57.90N  1 1.50E** 
61.55N  4.65E 
67.70N  33.70E 
61.  ION  30.20E 
59.30N  27.20E 
60.19N  5.25E 
58.90N  5.98E 
59.40N  28.50E 
60.60N  29.20E 
59.72N  22.56E 
59.40N  28.50E 
59.30N  27.60E 
59.40N  28.50E 
61. ION  30.20E 
61. ION  30.20E 
50.18N  12.35E 
50.17N  12.44E 
59.50N  25.00E 
60.90N  29.30E 
59.30N  27.60E 
60.0QN  28.50E** 
61.40N  31.60E 
59.40N  28.50E 
61. ION  30.20E 


EX-E7 

EQ 

EQ 

EQ 


Pn  Lg 
Pn  Lg 


Pn  Lg 
Pn  Lg 


EX-BLA 

Pn 

Lg 

1.9 

(R) 

EQ 

Pn 

l8 

2.8 

EX-V5 

Pn 

Lg 

2.8 

EX-BLA 

Pn 

Lg 

2.4 

EX-TIT 

Pn 

Lg 

2.1 

— 

Pn 

Lg 

2.0 

(R) 

EX-E3 

Lg 

2.6 

EX-E7 

Pn 

Lg 

2.7 

EX-TIT 

Pn 

Lg 

1.9 

EX-V12 

Pn 

Lg 

2.9 

EX-V8 

Pn 

2.9 

EX-E4 

Lg 

2.4 

(H) 

— 

Pn 

Lg 

1.9 

EX-E3 

Lg 

2.5 

EX-E8 

Lg 

3.0 

EQ 

Pn 

Lg 

3.0 

EX-V4 

Pn 

2.8 

HQ 

Pn 

Lg 

2.1 

(R) 

EQ 

Pn 

Lg 

3.0 

EX-K1 

Lg 

2.8 

EX-V2 

Pn 

<2.0 

(H) 

EX-E4 

Lg 

3.1 

— 

Pn 

Lg 

2.2 

EQ 

Pn 

Lg 

1.9 

EX-E8 

Lg 

3.2 

EX-VI  1 

Pn 

Lg 

2.0 

(R) 

EQ 

Pn 

Lg 

2.0 

EX-E8 

Pn 

Lg 

3.3 

EX-E6 

Lg 

2.8 

EX-E8 

Pn 

Lg 

2.8 

EX-V2 

Pn 

2.4 

(H) 

EX-V2 

Pn 

2.5 

(H) 

EQ 

Pn 

3.2 

(R) 

EQ 

Pn 

2.6 

(R) 

EX-E3 

Lg 

2.6 

EX-VIC 

Pn 

Lg 

2.9 

EX-E6 

Lg 

2.6 

EQ 

Pn 

Lg 

2.7 

EX-V4 

Pn 

2.2 

(H) 

EX-E8 

Pn 

Lg 

3.3 

EX-V2 

Pn 

2.4 

(H) 

r  Event 

Date 

Time 

Location 

Type 

Phases 

Ml 

$ 

5 

’•Lh 


§ 


46 

12-28-85 

47 

12-29-85 

48 

12-30-85 

49 

12-30-85 

50 

12-31-85 

51 

12-31-85 

52 

12-31-85 

53 

12-31-85 

54 

1-  3-86 

55 

1-  7-86 

56 

1-  7-86 

57 

1-  9-86 

58 

1-  9-86 

59 

1-13-86 

60 

1-15-86 

61 

1-16-86 

62 

1-17-86 

63 

1-17-86 

64 

1-20-86 

65 

1-21-86 

66 

1-25-86 

67 

1-25-86 

68 

1-31-86 

69 

1-31-86 

70 

1-31-86 

71 

1-31-86 

72 

2-  2-86 

73 

2-  3-86 

74 

2-  5-86 

75 

2-  5-86 

76 

2-  6-86 

77 

2-  6-86 

78 

2-  6-86 

79 

2  -7-86 

80 

2  -7-86 

81 

2-  7-86 

82 

2-  7-86 

83 

2  -7-86 

84 

2-10-86 

85 

2-13-86 

86 

2-14-86 

87 

2-14-86 

88 

2-15-86 

89 

2-16-86 

90 

2-16-86 

91 

2-17-86 

92 

2-18-86 

93 

2-18-86 

nK] 

K 


57.69N  26.54E 
73.29N  6.86E 
59.50N  26.50E 
59.31N  27.34E 
73.36N  6.77E 
73.29N  6.70E 
63.20N  27.80E 
58.34N  6.43E 
61.90N  30.60E 
60.92N  29.05E 
58.34N  6.43E 
66.80N  21.80E* 
59.30N  28.10E 
59.64N  24.07E 
59.40N  28.50E 
59.40N  28.50E 
59.30N  28.10E 
58.34N  6.43E 
50.18N  12.3  IE 
55.30N  13.60E* 
57. ION  7.00E** 
61.48N  16.94E 
65.39N  10.65E 
61.  ION  29.90E 
59.30N  28.10E 
58.34N  6.43E 
67.60N  34.00E 
59.60N  1.43E 
62.60N  6.80E* 
62.74N  4.50E 
62.90N  4.86E 
59.30N  28.10E 
67.  ION  20.60E 
64.70N  30.70E 
59.40N  28.40E 
59.20N  31.00E** 
67.60N  34.20E 
66.45N  14.89E 
59.40N  28.50E 
62.61N  5.07E 
58.34N  6.43E 
58.34N  6.43E 
59.86N  5.73E 
67.  ION  20.60E 
61.69N  4.90E 
59.30N  27.20E 
59.30N  27.20E 
64.70N  30.70E 


EQ 

EQ 

EX-E9 

EQ 

EQ 

EX-M7 

EX-TIT 

EX-V7 

EX-TIT 

EQ 

EX-E7 

EX-E8 

EX-E8 

EX-E7 

EX-TIT 

EQ 

EQ 

EQ 

EQ 

EQ 

EX-V8 

EX-E7 

EX-TIT 

EX-K2 

EQ 

EQ 

EQ 

EQ 

EX-E7 

EX-R1 

EX-VIO 

EX-E12 

EQ 

EX-K5 

EQ 

EX-E8 

EQ 

EX-TIT 

EX-TIT 

EQ 

EX-R1 

EQ 

EX-E4 

EX-E4 

EX-VIO 


Pn 

Lg 

2.8 

Pn 

4.7 

(P) 

Lg 

2.7 

Lg 

2.7 

Pn 

4.8 

(P) 

Pn 

4.6 

(P) 

Pn 

2.3 

(H) 

Pn 

Lg 

2.1 

Pn 

Lg 

2.9 

Pn 

Lg 

2.8 

Pn 

Lg 

1.8 

Lg 

2.6 

Pn 

Lg 

2.5 

Pn 

Lg 

2.7 

Pn 

Lg 

3.4 

Pn 

Lg 

2.7 

Pn 

Lg 

3.3 

Pn 

Lg 

2.3 

Pn 

Lg 

3.3 

Pn 

Lg 

2.5 

Pn 

Lg 

1.7 

Pn 

Lg 

2.9 

Pn 

Lg 

2.5 

Pn 

Lg 

3.3 

Pn 

Lg 

3.3 

Pn 

Lg 

1.9 

Lg 

2.9 

Lg 

1.9 

Lg 

1.6 

Pn 

Lg 

2.3 

Pn 

Lg 

1.9 

Pn 

Lg 

2.6 

Lg 

2.6 

Pn 

Lg 

3.1 

Lg 

2.5 

Lg 

1.9 

(R) 

Pn 

Lg 

2.8 

(H) 

Pn 

Lg 

2.2 

(H) 

Pn 

Lg 

2.7 

Lg 

2.2 

Pn 

Lg 

2.4 

Pn 

2.3 

(B) 

Lg 

1.8 

Lg 

2.5 

Pn 

Lg 

1.7 

Lg 

2.5 

(H) 

Lg 

3.1 

Pn 

Lg 

2.6 

Event 

Date 

Time 

Location 

94 

2-23-86 

6.14 

67.60N  34.00 E 

95 

2-26-86 

2:12 

62.76N  5.29E 

96 

3-  5-86 

12:13 

59.50N  26.50E 

97 

3-  5-86 

57.20N  7.00E** 

98 

3-  7-86 

59.30N  28.10E 

99 

3-  8-86 

16:21 

61.67N  2.58E 

62.81N  4.91E 

101 

3-10-86 

12:02 

59.30N  28. 10E 

102 

3-11-86 

12:02 

59.30N  28. 10E 

103 

3-13-86 

10:27 

61. ION  29.90E 

104 

3-13-86 

11:39 

60.70N  29.00E 

105 

3-14-86 

8:33 

67.60N  34.20E 

106 

3-21-86 

13:02 

59.50N  25.00E 

107 

3-24-86 

11:18 

59.30N  27.20E 

108 

3-25-86 

9:05 

62.76N  4.76E 

109 

3-27-86 

12:24 

59.40N  28.50E 

110 

3-30-86 

3:23 

61.66N  4.53E 

111 

4-  4-86 

13:13 

58.34N  6.43E 

112 

4-  4-86 

22:43 

70.86N  8.91E 

113 

4-  7-86 

0:35 

61.84N  4.88E 

114 

4-14-86 

14:55 

59.49N  24.1  IE 

115 

4-15-86 

10:53 

60.90N  29.30E 

116 

4-16-86 

11:51 

60.39N  5.34E 

117 

4-16-86 

13:15 

58.15N  5.97E 

118 

4-18-86 

8:33 

67.60N  34.00E 

119 

4-19-86 

10:59 

61.  ION  30.20E 

120 

4-28-86 

15:53 

60.18N  4.88E 

121 

4-29-86 

17:48 

59.82N  24.06E 

122 

4-30-86 

10:19 

59.31N  6.95E 

123 

5-  8-86 

17:14 

58.70N  17.99E 

124 

5-16-86 

15:02 

61.90N  30.60E 

125 

5-17-86 

16:01 

62.94N  4.94E 

126 

5-21-86 

8:57 

61.65N  31.38E 

127 

5-27-86 

18:36 

59.31N  6.95E 

128 

5-28-86 

17:52 

59.31N  6.95E 

129 

6-  3-86 

11:04 

59.17N  5.66E 

130 

6-  3-86 

14:30 

61.46N  4.08E 

131 

6-  4-86 

9:07 

61.50N  30.40E 

132 

6-  6-86 

13:14 

58.34N  6.43E 

133 

6-  7-86 

12:13 

59.20N  27.60E 

134 

6-12-86 

9:31 

61.50N  30.40E 

135 

6-12-86 

13:04 

60.90N  29.30E 

136 

6-13-86 

14:41 

59.66N  24.28E 

137 

6-15-86 

15:01 

61.67N  3.85E 

138 

6-16-86 

15:59 

60.04N  7.24E 

139 

6-19-86 

3:55 

59.3 IN  6.95E 

140 

6-20-86 

22:08 

61.47N  3.92E 

141 

6-23-86 

13:14 

58.34N  6.43E 

I 


Type  Phases 


EX-K2 

Pn 

Lg 

EQ 

Pn 

Lg 

EX-E9 

Lg 

EQ 

Pn 

Lg 

EX-E7 

Pn 

Lg 

EQ 

Pn 

Lg 

EQ 

Pn 

Lg 

EX-E7 

Pn 

Lg 

EX-E7 

Pn 

Lg 

EX-V8 

Pn 

EX-V5 

Pn 

Lg 

EX-K5 

Lg 

EX-E3 

Lg 

EX-E4 

Pn 

Lg 

EQ 

Pn 

Lg 

EX-E8 

Pn 

Lg 

EQ 

Pn 

Lg 

EX-TIT 

Pn 

Lg 

EQ 

Pn 

EQ 

Pn 

Lg 

-- 

Pn 

Lg 

EX-VIC 

Pn 

Lg 

EX-NYG 

Lg 

- 

Pn 

Lg 

EX-K2 

Lg 

EX-V2 

Pn 

Lg 

- 

Pn 

Lg 

EQ 

Lg 

EX-BLA 

Pn 
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Figure  4.3.  Locations  of  mines  listed  in  Table  4.2. 


SAIAP  also  uses  a  broad  band,  frequency-wavenumber  calculation  (Kvaerna  and  Ringdal,  1986) 
that  gives  more  accurate  azimuth  estimates  than  those  used  in  the  RONAPP  locations. 

The  Mi  are  based  on  the  Lg  amplitude  computed  by  RONAPP.  They  differ  slightly  from  the 
RONAPP  Mi  in  the  NORESS  bulletin  by  being  distance-conected  to  the  event  location  com¬ 
puted  by  one  of  the  independent  networks.  For  some  events  RONAPP  had  no  deiected  Lg 
phase,  or  chose  the  wrong  phase  as  Lg.  For  these  events  the  Mi  is  that  reported  by  one  of  the 
independent  networks  (if  available)  or  the  RONAPP  uncorrected  magnitude.  These  are  listed  in 
Table  4.1  for  information  but  were  not  included  in  any  subsequent  analysis.  The  relation 
between  seismic  moment  and  magnitude  (Section  5.4)  was  derived  using  only  events  with 
distance-corrected  NORESS  magnitudes. 

The  190  events  used  in  the  inversion  include  109  explosions,  65  presumed  earthquakes,  and  16 
events  of  unknown  source  type.  The  Pn  inversion  included  152  of  these  events,  with  83  explo¬ 
sions,  56  presumed  earthquakes,  and  13  events  of  unknown  source  type.  The  Lg  inversion 
included  170  events,  including  97  explosions,  58  presumed  earthquakes,  and  15  unknown.  Of 
the  190  events,  132  were  used  for  both  the  Pn  and  the  Lg  analyses.  The  NORESS  array  loca¬ 
tion  and  epicenters  of  events  used  in  the  inversion  are  plotted  in  Figure  4.4.  The  distance 
range  300-450  km  is  dominated  by  events  in  west  to  southwest  Norway,  while  events  from 
700  to  1200  km  are  located  primarily  to  the  east  of  NORESS.  If  the  attenuation  along  paths  to 
the  east  is  distinctly  different  from  that  along  paths  to  the  west,  it  would  not  be  possible  to 
combine  them  in  the  same  inversion.  However,  separate  inversions  were  run  for  restricted 
azimuth  windows,  and  significant  path  differences  were  not  observed.  From  this  we  conclude 
that  all  of  the  data  are  adequately  represented  by  a  single  frequency-dependent  Q  model. 

Appendix  A  lists  an  additional  68  events  that  were  processed  but  that  were  not  used  in  the 
inversion.  Most  were  excluded  because  of  low  signal/noise  ratios.  However,  there  is  a  set  of 
southwest  Norway  off-shore  events  that  were  excluded  because  their  spectra  were  inconsistent 
with  other  events  at  approximately  the  same  distance  (Henson  and  Bache,  1986).  It  is  not  clear 
whether  the  difference  is  due  to  anomalous  path  effects  or  to  a  different  source  spectrum. 
These  events  are  discussed  in  more  detail  in  Section  6.5. 


Figure  4.4.  Events  used  in  the  generalized  inversion  for  (a)  Pn  and  (b)  Lg.  Earthquakes, 
explosions,  and  sources  of  unknown  type  are  plotted  with  different  symbols.  The 
Lg  inversion  used  spectra  of  170  events  and  the  Pn  inversion  included  152 
events. 


5.  ATTENUATION  AND  SEISMIC  MOMENT  ESTIMATES 


In  this  section  we  present  the  primary  results  of  this  study.  The  generalized  inversion 
described  in  Section  3  is  applied  to  regional  Pn  and  Lg  spectra  recorded  at  NORESS  from  190 
events  with  epicentral  distances  ranging  between  200  and  1300  km.  The  data  used  for  the 
separate  Pn  and  Lg  inversions  were  described  in  Section  4.  In  this  section,  we  adopt  the 
Mueller-Murphy  explosion  source  model  (3.1.4)  and  geometric  spreading  defined  by  (3.1.2) 
with  m  -  1/2  and  rg  *=  100  km  for  Lg  and  m  =  1.3  and  rg  =  1  km  for  Pn.  In  Section  6,  the 
dependence  of  our  results  on  these  assumptions  is  investigated. 

For  a  given  set  of  source  and  spreading  assumptions,  the  inversion  clearly  defines  a  broad 
minimum  in  the  data  residuals  corresponding  to  a  suite  of  models  that  fit  the  data  equally  well 
in  a  least-squares  sense.  These  models  involve  trade-offs  among  Mq,  Qq,  and  T|.  The  low  fre¬ 
quency  spectral  level  depends  upon  Mq  and  Qq.  If  Mq  increases,  then  Qq  will  decrease  to 
preserve  the  fit  to  the  long  period  spectral  level.  This  trade-off  is,  however,  limited  by  the 
spectral  slope  at  low  frequencies.  In  response  to  the  high  frequency  spectrum,  decrease  in  Qq 
trades-off  with  increase  in  T).  The  parameter  c  relating  Mq  to  corner  frequency  (3.1.6)  is  fixed 
by  the  few  events  large  enough  to  have  a  comer  frequency  within  the  band  used  in  the  inver¬ 
sion.  Therefore,  when  Mq  increases,  c  increases  to  retain  approximately  the  same  corner  fre¬ 
quency  for  the  larger  events.  To  resolve  these  trade-offs  among  models  that  give  essentially  the 
same  data  variance,  we  add  the  constraint  that  the  derived  source  parameters  for  Lg  and  Pn, 
which  are  inverted  separately,  be  consistent.  In  this  section,  the  results  for  our  "preferred 
model"  are  presented.  The  trade-off  among  model  parameters  is  discussed  in  detail  in  Section 
6. 

5.1  Regional  Lg  spectra 

Regional  Lg  spectra  of  170  events  (Table  4.1)  were  inverted  simultaneously  for  seismic 
moment  and  Q(f).  At  ranges  greater  than  about  800  km,  the  Lg  signal/noise  ratio  is  inadequate 
for  most  events  above  7  or  8  Hz,  while  at  shorter  distances  Lg  spectra  are  probably  contam¬ 
inated  by  Sn  coda  at  high  frequency  (Chun,  et  al.,  1987;  Ringdal,  1986;  Shin  and  Herrmann, 
1987).  The  Lg  inversion  was  therefore  band  limited  to  frequencies  between  1  and  7  Hz.  The 
spectra  were  smoothed  over  a  2  Hz  bandwidth  and  sampled  every  0.25  Hz.  The  number  of 
data  in  the  Lg  inversion  is  then  4250,  and  the  number  of  parameters  is  173. 


The  Lg  Q  resulting  from  the  inversion  is  Q(f)  =  350/ °'41.  Our  short  time  windows  include 
the  onset  of  Lg  which  consists  of  higher  mode  surface  waves  sampling,  predominantly,  the 
lower  crust.  Our  QLg(f)  is  therefore  an  approximation  to  Q  of  the  lower  crust,  although  it  also 
includes  the  effects  of  apparent  attenuation  due  to  scattering.  The  parameter  c  derived  from  Lg 
spectra  is  28.7.  The  inverted  long  period  source  levels  for  each  event  are  tabulated  in  Appendix 
B.  Selected  examples  of  the  fit  of  theoretical  spectra  to  data  spectra  are  shown  in  Figure  5.1. 
All  of  these  events  have  comer  frequencies  greater  than  6  Hz,  so  the  source  parameterization 
has  only  a  minor  effect  on  the  results.  A  complete  catalog  of  the  fit  to  all  spectra  used  in  the 
inversion  can  be  found  in  Appendix  B. 

Table  2.1  and  Figure  2.1  summarize  reported  Lg  attenuation  functions  for  various  regions.  It 
was  noted  in  Section  2  that  eastern  North  America  (ENA)  is  characterized  by  a  higher  Q0  and 
lower  frequency  exponent  than  the  western  Untied  States  (WUS).  Our  Qig{f)  is  between  the 
ENA  and  WUS  estimates.  This  is  surprising  since  one  would  expect  attenuation  in  Scandina¬ 
via  would  be  more  similar  to  the  tectonically  stable  North  American  shield.  However, 
predicted  spectral  shapes  for  our  Q  model  are  more  similar  to  those  for  the  ENA  estimates 
than  those  based  on  WUS  Q  values.  This  is  illustrated  in  Figure  5.2  which  plots 
exp (~nftlQ(f))  for  the  Q  functions  listed  in  Table  2.1,  and  for  our  Q(J).  In  particular,  the 
spectral  shape  between  1  and  7  Hz  for  our  Q  model  is  similar  to  that  predicted  for  the  low 
exponent  ENA  Q  models  (entries  6  and  8  of  Table  2.1).  On  the  other  hand,  a  large  Q0  and 
strong  frequency-dependence  (entry  10)  result  in  spectra  that  severely  overestimate  the 
observed  ratio  of  high  to  low  frequency  spectral  content  of  NORESS  data.  Our  (2z.g(/)  is  actu- 
ally  very  similar  to  that  representing  central  France  (Campillo,  et  al.,  1985).  The  vertical  offset 
between  the  ENA  curves  and  ours  indicates  that  our  model  predicts  considerably  lower  ampli¬ 
tudes  for  events  of  equivalent  Mq.  It  is  not  clear  whether  paths  to  NORESS  are  distinctly 
different  from  ENA  paths,  or  if  the  disparate  Q  values  can  be  attributed  to  different  methodolo¬ 
gies. 

5.2  Regional  Pn  spectra 

Regional  Pn  spectra  of  152  events  (Table  4.1)  were  inverted  between  1  and  15  Hz,  smoothed 
over  a  2  Hz  bandwidth,  and  sampled  every  0.25  Hz.  However,  we  note  that  at  ranges  less 
than  300-400  km,  Pn  and  Pg  arrive  within  the  same  5s  window  and  it  is  possible  that  some  of 
our  short  range  spectra  include  a  Pg  contribution.  The  number  of  data  in  the  Pn  inversion  is 
8664  and  the  number  of  parameters  is  155.  The  Pn  spectra  were  corrected  for  geometric 
spreading  and  inverted  for  seismic  moment  and  Q(J). 


Frequency  (Hz) 


Figure  5.2.  Lg  spectral  attenuation  for  the  Q  models  listed  in  Table  2.1  and  for  our  Scandi¬ 
navian  Q  estimate.  The  curves  are  identified  by  their  numbered  entries  in  Table 
2.1.  Attenuation  is  plotted  as  exp(-7t/f/Q(/))  where  the  travel  time  t  is  appropri¬ 
ate  for  a  range  of  1000  km.  Note  the  distinct  offset  between  eastern  North 
America  (ENA)  and  western  United  States  (WUS). 


The  apparent  attenuation  from  the  Pn  inversion  is  represented  by  Q(f)  =  300/ 049.  The 
comer  frequency  scaling  parameter  c  derived  from  Pn  spectra  is  29.2.  The  inverted  long  period 
source  levels  for  each  event  are  tabulated  in  Appendix  B.  Selected  examples  of  the  fit  of 
theoretical  spectra  to  data  spectra  are  shown  in  Figure  5.3.  The  Pn  amplitude  is  more  variable 
than  the  Lg  amplitude  probably  due  to  greater  sensitivity  to  source  radiation  pattern, 
focussing/defocussing,  and  scattering.  Nevertheless,  acceptable  fits  to  the  Pn  spectra  were 
achieved  with  our  simple  parameterization.  A  complete  catalog  of  the  fits  to  all  the  spectra 
used  in  the  inversion  can  be  found  in  Appendix  B. 


5  J  Pn/Lg  consistency 

An  important  constraint  on  the  attenuation  estimates  is  the  consistency  of  the  derived  source 
parameters  from  the  independently  inverted  Pn  and  Lg  spectra.  In  particular,  the  ratio  of  the 
Lg  to  Pn  long  period  source  levels  should  be  range-independent  Range  dependence  would 
indicate  that  geometric  spreading  and/or  Q  have  been  improperly  modeled  for  one  or  both 
phases.  Figure  5.4  is  a  plot  of  So(Lg)/SQ(Pn)  for  the  132  events  common  to  both  inversions. 
Note  that  for  explosions  the  ratio  does  not  show  evidence  of  range  dependence.  Examples  of 
models  that  demonstrate  an  obvious  range  dependence  are  presented  in  Section  6.  There  is  a 
much  larger  scatter  in  the  earthquake  ratios,  as  expected  from  source  radiation  pattern  effects. 
Based  on  the  explosions,  the  range-independent  value  of  SovC^/Sq^h)  is  approximately  0.7. 
Combining  equations  (3.1.7),  (3.1.9)  and  (3.1.10),  this  implies  for  earthquakes  and  explosions 
of  equal  moment 


P  fit 
PcPc 


0.7 


(5.3.1) 


Assuming  surface  values  of  p,  -  2.5  gm/cm3  and  Cts  =  5.0  km/s  and  average  crustal  values  of 
pc  -  2.7  gm/cm3  and  (3C  «  3.5  km/s,  we  have  1^0.26.  Simply  stated,  this  means  that  the 
average  Lg  earthquake  excitation  is  approximately  4  times  the  average  Lg  explosion  excitation 
for  sources  of  equal  moment.  We  note  that  this  is  consistent  with  earlier  studies.  For  example, 
Willis  (1963),  compared  Lg  amplitudes  from  an  NTS  explosion  and  a  co-Iocated  earthquake  at 
a  range  of  450  km  and  found  a  factor  of  5  for  their  ratio.  Pomeroy  (1977)  found  a  ratio  of  3-5 
for  relative  Lg  excitation  by  comparing  12  earthquakes  to  the  SALMON  nuclear  explosion 
detonated  in  Mississippi.  However,  other  studies  have  found  a  less  distinct  separation  (e.g., 
Murphy  and  Bennett,  1982;  Nuttli,  1981).  We  will  address  the  sensitivity  of  K  with  respect  to 
spreading  assumptions  in  Section  6.3.  Finally,  we  note  in  Figure  5.4  that  the  50(Z.^)/50(/>«) 
ratio  separates  many  of  the  earthquakes  from  the  explosions,  however,  there  is  considerable 
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Figure  3.3.  Selected  comparisons  of  theoretical  Pit  spectra,  based  on  the  inversion  results,  to 
observed  spectra.  Pn  spectra  of  earthquakes  (events  137,  64  and  47  of  Table 
4.1)  are  on  the  left  and  Pn  spectra  of  mine  blasts  (events  102,  30,  and  82  of 
Table  4.1)  are  on  the  right  The  explosion  spectra  are  labeled  by  the  Helsinki 
mine  identification  code  (Table  4.2). 
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Figure  5.4.  Ratio  of  Lg  to  Pn  inverted  long  period  source  levels  as  a  function  of  epicentral 
distance.  Explosions,  earthquakes  and  events  of  unknown  source  type  are  plotted 
with  different  symbols. 
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overlap  between  the  two  populations.  The  large  ratios  for  some  earthquakes  are  presumably 
due  to  the  source  radiation  pattern  for  Pn. 


Corner  frequency  is  not  well-constrained  in  the  Lg  inversion.  Most  of  the  events  used  in  this 
study  have  local  magnitudes  less  than  3.0  and  correspondingly  high  corner  frequencies.  The  1- 
7  Hz  bandwidth  used  in  the  Lg  inversion  was  insufficient  to  resolve  the  corner  frequency  with 
much  confidence.  We  compute  the  ratio  of  the  Pn  derived  corner  frequency  to  the  Lg  derived 
corner  frequency  from 


_  c(Pn)  S0(Lg) 
fc(Lg)  c(Lg)  S0(Pn) 


With  c(Pn)  -  29.2,  c(Lg)  -  28.7  and  Sp(Lg)/  Sp(Pn)  ratio  of  0.7  typical  for  the  explosions, 
the  corner  frequency  ratio  is  0.9,  which  is  not  significantly  different  from  unity.  We  conclude 
that  the  Q(f)  models  we  have  derived  result  in  consistent  source  parameters  for  the  events 
common  to  both  inversions,  which  substantially  improves  confidence  in  the  validity  of  the 
results. 


5.4  Seismic  moments 

Seismic  moment  is  estimated  from  the  inversion  parameter,  S0,  using  equations  (3.1.7)  to 
(3.1.10).  The  events  are  at  different  depths  and  locations,  so  our  estimates  are  "relative 
moments"  in  that  we  use  the  same  near-source  material  properties  for  all  events.  The  near¬ 
surface  values  used  are  ps  *  2.5  gm/cm3  and  (xs  -  5.0  km/s,  and  crustal  values  used  are  pc  = 
2.7  gm/cm3  and  Pc  -  3.5  km/s.  Appendix  B  tabulates  the  estimated  moments  for  most  of  the 

events  used  in  the  inversion.  Moments  were  not  estimated  for  earthquakes  for  which  only  the 

Pn  data  were  inverted,  or  for  events  of  unknown  source  type.  Similarly,  moments  were  not 
estimated  for  explosions  for  which  only  Lg  data  were  inverted  because  of  uncertainty  in  the 
constant,  K  (3.1.10). 

In  Figure  5.5  we  plot  explosion  moment  versus  NORESS  local  magnitude.  The  least-squares 
linear  fits  to  LogM0  are  given  by 

Pn  :  LogMpxp  =1.12  +  17.6  (5.4.1) 

Lg  :  LogM'p'P  =  1.16  ML  +  17.5  (5.4.2) 


Inverted  explosion  moment  versus  NORESS  local  magnitude  for  (a)  Pn  (3.1.7) 
and  (b)  Lg  (3.1.10)  with  K  -  0.26.  The  best  fitting  straight  line  to  LogM$  is 
shown  in  each  figure. 


For  Lg,  we  assume  K  -  0.26.  The  tc  effects  the  intercept,  but  not  the  slope  in  (5.4.2).  Our 
distance-corrected  NORESS  M i  were  used  because  they  provide  a  consistent  magnitude  meas¬ 
ure  for  the  data  set.  Events  for  which  we  do  not  have  this  M ^  were  not  included.  Comparing 
the  two  moment-magnitude  relationships,  the  Mqxp  derived  from  Lg  is  equal  to  the  Mqxp  from 
Pn  when  ML  •  2.5.  This  is  near  the  center  of  the  ML  distribution  for  our  data  set  and  pro¬ 
vides  further  confirmation  of  the  internal  consistency  of  our  analysis. 

The  earthquake  moments  were  estimated  from  Lg  spectra  using  (3.1.9)  and  the  result  is 
displayed  in  Figure  5.6.  The  least-squares  linear  fit  to  these  data  is 

LogMtf  =  1.04  Ml  +  17.1  (5.4  3) 

Comparing  (5.4.3)  to  (5.4.1),  for  equivalent  moments  the  earthquake  M \  (i.e.,  log  Lg  ampli¬ 
tude)  is  larger  than  the  explosion  by  0.6  to  0.7,  or  a  factor  of  4  to  5  in  Lg  amplitude. 
Table  5.1  compares  our  earthquake  moment-magnitude  relation  to  results  of  previous  studies. 
The  Bungum,  et  al.  (1982)  study  used  near-field  S  wave  spectra  from  the  1978  Meloy  earth¬ 
quake  sequence  in  northern  Norway.  In  general,  our  results  for  earthquake  moment  as  a  func¬ 
tion  of  local  magnitude  are  in  agreement  with  those  of  similar  studies  and,  in  particular,  with 
the  result  of  Bungum,  et  al.  (1982)  which  utilized  data  horn  the  same  geographic  region.  In 
that  magnitude  does  not  enter  the  inversion,  this  result  lends  considerable  support  to  our 
derived  Q  model. 

5.5  Corner  frequency 

Cube-root  corner  frequency  scaling  (3.1.6)  was  assumed  and  we  inverted  for  the  constant,  c, 
relating  corner  frequency  to  long  period  source  level.  Appendix  B  tabulates  corner  frequencies 
for  each  event  using  the  c(Pn)  and  c(Lg)  from  the  inversion  and  the  Sq  for  each  event.  On 
average,  our  results  indicate  that  an  M i  =  3.0  earthquake  has  a  corner  frequency  of  approxi¬ 
mately  1 1  Hz.  Most  of  the  events  have  M i  <  3.0  and  thus  have  high  corner  frequencies. 
Therefore,  the  data  do  not  clearly  resolve  source  corner  frequency,  but  we  note  that  our  results 
are  consistent  with  those  from  similar  studies.  For  example,  Bungum,  et  al.  (1982)  expressed 
the  corner  frequency-moment  relationship  for  the  1978  Meloy  earthquake  sequence  in  northern 
Norway  as 


fc  =  -205  +  30.8 LogM0  -  ( Log\10 )2  (5.5.1) 
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Table  5.1.  Moment-magnitude  relations. 


-0.4  - 2.2 

1.0  -  3  5 

2.1  -  4.3 
10  -  4.2 
18  -  4.9 

1.1  -  3  8 


p  (gm/cm  ) 

p  (km/s) 

LogM0 

Distance 

2.7 

3.5 

0.90Ml  +  17.5 

near-field 

2.7 

3.5 

l.OlA/i  +  16. 7t 

near-field 

2.8 

3.8 

l.l8Aft  +  16. 6t 

regional 

2.8 

3.8 

0.94ML  +  17.32 

regional 

2.5 

3.5 

1.26 Ml  +  16.44 

<  200  km 

2.7 

3.5 

\.0AMl  +  17.1 

regional 

_ Reference 

Bungum,  et  al.  (1982) 

Mueller  and  Cranswick  (1985) 
Shin  and  Herrmann  (1987) 
Hasegawa  (1983) 

Dwyer,  at  al.  (1983) 

Sereno,  et  al.  (this  study) 


+  We  estimated  the  regression  coefficients  from  tabulated  moments  and  magnitudes. 


Adopting  their  moment-magnitude  relation  (Table  5.1),  this  gives  a  corner  frequency  of  9  Hz 
for  an  M L  -  3.0  event.  Similarly,  assuming  cube-root  scaling,  a  least-squares  fit  to  the  Mirami- 
chi  aftershock  comer  frequencies  derived  by  Mueller  and  Cranswick  (1985)  gives  a  corner  fre¬ 
quency  of  12  Hz  for  an  ML  -  3.0  earthquake.  In  general,  the  comer  frequencies  derived  in  this 
study  are  within  a  few  Hertz  of  the  near-field  estimates. 


6.  ERROR  ANALYSIS 


The  separate  Pn  and  Lg  inversions  are  each  characterized  by  a  broad  minimum  in  the  data 
residuals,  and  therefore  a  range  of  models  that  fit  the  data  equally  well.  In  Section  S  we 
described  our  "preferred  model"  and  in  this  section  we  define  the  range  of  acceptable  models 
and  the  basis  upon  which  our  selection  was  made.  The  following  list  summarizes  the  selection 
process. 

(1)  Identify  the  full  range  of  parameter  trade-offs  for  each  phase.  For  fixed  source  and 
spreading  assumptions,  the  separate  Pn  and  Lg  inversions  identify  the  range  of  Q 
models  that  produce  essentially  the  same  data  variance. 

(2)  Select  Pn  and  Lg  models  that  give  consistent  source  parameters  for  the  events  com¬ 
mon  to  both  inversions.  That  is,  constrain  So(Lg)tSo(Pn)  to  be  range-independent 
This  reduces  the  acceptable  solutions  to  a  set  of  model  pairs.  That  is,  for  a  fixed 
Lg  Q  model,  tight  bounds  are  placed  on  acceptable  Pn  Q  models.  From  the  range- 
independent  value  of  the  source  level  ratio,  estimate  the  relative  Lg  excitation  of 
earthquakes  and  explosions. 

(3)  Analyze  the  trade-off  between  data  variance  and  Mq-Ml  variance.  It  was  found  that 
models  which  reduced  the  data  variance  increased  the  variance  of  the  Mq-M ^  rela¬ 
tion.  By  excluding  models  that  increase  either  variance  without  significant  decrease 
in  the  other,  we  define  the  range  of  acceptable  model  pairs.  The  "preferred  model" 
(for  fixed  source  and  spreading  assumptions)  is  selected  from  the  mid-range  of  the 
acceptable  model  pairs. 

(4)  Change  source  assumptions.  It  was  found  that  our  results  were  not  strongly  depen¬ 
dent  upon  the  details  of  our  source  assumptions  because  of  the  high  corner  frequen¬ 
cies  for  most  of  the  events. 

(5)  Change  spreading  assumptions  We  fix  Lg  spreading  (which  is  well-constrained 
theoretically  and  empirically)  and  change  the  Pn  spreading  rate.  The  new  set  of  Pn 
Q  models  are  input  to  steps  (2)  and  (3)  above  to  define  acceptable  models  for  this 
Pn  spreading  rate  Next,  compare  the  Lg  earthquake-explosion  ratio  to  the  results  of 
other  studies  to  bound  the  Pn  spreading  rate. 


We  present  in  Section  6.1  the  formal  inversion  errors  in  the  form  of  a  parameter  covariance 
matrix.  However,  we  note  that  the  relatively  small  formal  errors  obtained  do  not  reflect  the 
true  uncertainty  in  our  parameter  estimates.  In  Section  6.2  we  discuss  the  trade-offs  among 
model  parameters  for  the  fixed  spreading  and  source  assumptions  used  in  Section  5.  The 
consequences  of  changing  these  assumptions  are  explored  in  Section  6.3.  Section  6.4  describes 
observed  systematic  variations  in  source  moment  versus  Lg  amplitude  for  different  mines. 
Section  6.5  presents  examples  of  data  that  were  not  well-modeled  by  the  inversion  results. 

6.1  Formal  inversion  errors 

The  variance  of  the  inversion  is  defined  as  the  sum  of  the  squared  data  residuals  divided  by  the 
number  of  data  (nd)  minus  the  number  of  parameters  ( np ).  That  is 

Oj  = ’  IW,  -  ^)2  (6-1.1) 

nd-np  t=, 

In  this  case,  the  data  are  log  amplitudes  corrected  for  geometric  spreading.  Figure  6.1  displays 
the  Pn  and  Lg  spectra  for  a  magnitude  3.3  German  earthquake  at  a  range  of  1 170  km.  Super¬ 
imposed  are  theoretical  spectra  for  the  model  presented  in  Section  5  bounded  by  single  stan¬ 
dard  deviation  curves  (Od).  The  variance  of  the  Pn  inversion  is  approximately  0.06  and  the  Lg 
variance  is  approximately  0.03.  Note  that  Oj  is  the  variance  of  the  smoothed  spectra  and  is 
therefore  lower  than  the  true  data  variance. 


Formal  error  estimates  on  the  individual  inversion  parameters  are  obtained  from  the  parameter 
covariance  matrix.  For  the  least-squares  solution,  (3.1.24),  the  parameter  covariance  matrix  is 


(AtA)-' 


(6.1.2) 


The  parameter  covariance  matrix  for  the  damped  least-squares  solution,  (3.1.25),  is 


(ATA  +  *  D~XATA(AtA  +  \I)~l 


(6.1.3) 


The  Lg  Q  model  presented  in  Section  5  had  Q$  *  349  with  a  formal  inversion  error  of  ±  8 
and  T\  =  0.408  ±  0.008.  The  Pn  Q  model  had  Q0  =  301  ±  9  and  r|  =  0.486  ±  0.0 10.  The 
corner  frequency  parameters  are  c(Pn )  =  29.24  ±  0.04  and  c(Lg)  =  28.66  ±  0.05.  The 
artificially  small  formal  errors  on  c  are  a  consequence  of  the  poor  resolution  of  this  parameter. 
That  is,  damping  has  decreased  the  variance  of  the  parameter  estimate  with  an  associated 
decrease  in  its  resolution.  The  variance  of  the  individual  Sq  estimates  increase  with  increasing 
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epicentral  distance  but,  on  the  average,  the  S0  standard  deviations  were  of  the  order  of  10%. 
The  standard  errors  of  the  explosion  moment-magnitude  relations  are  LogM 0  = 
112(  ±  0.08 )Ml  +  17. 6(  ±  0.2)  for  Pn  and  LogM0  =  1.16(  ±  0.08)ML  +  17.5(  ±  0.2)  for 
Lg.  The  Lg  earthquake  moment-magnitude  relation  is  LogM0  =  1.04(  ±  0.04)ML  + 
17.1(±0.1).  These  errors  must  be  taken  into  account  in  analyzing  the  predictability  of 
regional  wave  spectra  of  events  of  a  given  magnitude  (Section  7).  The  single  standard  devia¬ 
tion  curves  in  Figure  6.1  reflect  only  the  estimated  uncertainty  in  the  spectral  estimate  for  fixed 
moment,  and  do  not  include  the  scatter  in  the  moment-magnitude  relationship. 


6.2  Parameter  trade-off 

In  this  section,  we  discuss  the  tradeoffs  among  M0,  Q0,  T},  and  c  for  fixed  source  and  spread¬ 
ing  assumptions.  The  effect  of  varying  the  source  and  spreading  assumptions  will  be  discussed 
in  the  next  section.  The  low  frequency  spectral  level  depends  on  M0  and  Q0.  When  M0 
increases,  Qq  decreases  within  a  range  of  models  that  preserve  the  fit  to  the  long  period  spec¬ 
tral  level  and  slope.  In  response  to  the  high  frequency  spectrum,  when  Q0  decreases,  r[ 
increases.  The  corner  frequency  scaling  parameter  c  is  constrained  by  a  few  large  events,  and 
increases  in  the  estimates  for  M0  cause  an  increase  in  c  to  retain  approximately  the  same 
corner  frequency  for  these  events.  These  trade-offs  are  illustrated  in  Table  6.1  with  6  models 
( Qo  varying  between  300  and  427)  obtained  from  the  Lg  inversion  that  have  data  variances 
that  differ  by  less  than  5%. 

The  trade-offs  in  the  Pn  results  are  illustrated  in  Table  6.2  with  7  Pn  models  that  have  data 
variances  within  13%.  All  were  done  with  the  r~1,3  spreading  assumption.  If  the  Q  model  for 
Lg  is  fixed  then  the  requirement  that  SQ(Lg)!So(.Pn)  be  range-independent  places  tight  bounds 
on  acceptable  Pn  Q  models.  For  example,  in  Figure  6.2  we  plot  this  ratio  for  Lg  Qq  -  350 
and  three  Pn  models  with  Qq  -  401,  300,  and  201.  Only  explosions  were  included  in  this 
figure  because  the  earthquake  ratio  is  contaminated  by  unknown  Pn  radiation  pattern  effects. 
The  ratios  for  the  Qq  «  401  and  201  models  display  a  clear  range  dependence.  Thus,  if  Lg  Qq 
is  350,  then  the  corresponding  choice  for  Pn  Q0  is  300  and  these  are  the  models  used  in  Sec¬ 
tion  5. 

For  other  choices  of  Lg  Qq,  other  Pn  Q  models  are  consistent  with  a  range- independent  S0 
ratio.  In  particular,  Lg  Qq  =  400  suggests  Pn  Qq  ~  375  and  Lg  Q0  =  300  suggests  Pn 
Qq  ~  250.  There  is  no  obvious  basis  for  choosing  among  the  model  pairs.  However,  there  is 
a  trade-off  between  data  variance  and  variance  of  the  moment-magnitude  relation.  Figure  6.3 
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Table  6.2.  Pn  parameter  trade-off 


Model  Index 

Qo 

mm 

c 

Lo&M?  p 

Vd 

1 

498 

0.39 

16.0 

0.99Ml  +  17.63 

0.0703 

2 

401 

0.43 

20.8 

1.04A/l  +  17.63 

0.0667 

3 

373 

0.44 

22.8 

1.05Ml+  17.62 

0.0657 

4 

356 

0.45 

23.3 

l.01ML  +  17.61 

0.0653 

5t 

300 

0.49 

29.2 

1.12Ml  +  17.57 

0.0637 

6 

248 

0.53 

36.6 

1 . 1  9Ml  +  17.52 

0.0625 

7 

201 

0.57 

48.8 

1.29  Ml  +  17.43 

0.0614 

t  The  model  discussed  in  Section  5. 
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Figure  6.2.  Ratio  of  Lg  to  Pn  inverted  long  period  explosion  source  levels  as  a  function  of 
epicentral  distance.  Lg  Qq  is  fixed  for  all  3  plots  at  350  (#4  of  Table  6.1).  Pn 
Qn  is  varied  from  401  to  201  (numbers  2,  5,  and  7  of  Table  6.2). 


Trade-off  between  data  variance  (CT^)  and  variance  of  the  moment-magnitude 
relation  for  (a)  Pn  and  (b)  Lg.  Symbols  define  models  for  different  iterations  of 
the  inversion.  Three  of  the  models  are  identified  in  each  plot  by  their  Q0  values. 
For  Lg  these  are  models  2,  4,  and  6  of  Table  6.1  and  for  Pn  they  are  models  3, 
5,  and  6  of  Table  6.2. 


plots  this  trade-off  for  the  range  of  models  defined  in  Tables  6.1  and  6.2.  The  figure  demon¬ 
strates  that  the  lower  Q  models  reduce  the  data  variance  but  result  in  increased  scatter  in  the 
moment- magnitude  relationship.  On  the  other  hand,  the  higher  Q  models  reduce  the  Mq-Ml 
variance  but  increase  the  data  variance. 

From  Figure  6.3,  we  note  that  models  with  Pn  Qq  >  375  cause  a  rapid  increase  in  data  vari¬ 
ance  with  only  a  slight  reduction  in  Mq-M ^  variance.  Also,  models  with  Pn  Qq  <  250  increase 
the  Mq-Ml  variance  without  reduction  in  data  variance.  We  therefore  define  the  range  of 
acceptable  Pn  Q  models  as  Qo(Pn)  =  250-375  with  T]  between  0.53  and  0.44.  Therefore,  the 
range-independence  of  the  source  level  ratios  requires  the  range  of  acceptable  Lg  Q  models  to 
be  Qo(Lg)  =  300-400  with  T]  between  0.46  and  0.37  These,  of  course,  cannot  be  chosen  in  any 
combination.  For  a  fixed  Lg  Q  model,  a  tight  bound  on  Pn  Q  is  placed  by  the  range- 
independence  of  the  Lg  to  Pn  source  level  ratio.  From  the  mid-range  of  acceptable  models, 
we  selected  the  Lg  Qq  =350  and  Pn  Qq  =  300  model  pair  as  our  "preferred  model”  to  present 
in  Section  5. 

6.3  Result  dependence  on  parameterization 

In  Section  6.2,  we  discussed  the  parameter  trade-off  for  a  fixed  set  of  source  and  geometric 
spreading  assumptions.  In  this  section,  we  analyze  the  dependence  of  our  results  on  these 
assumptions. 

6.3.1  Source  parameterization 

All  of  the  results  of  this  study  were  obtained  by  assuming  an  to2  source  function  with  cube- 
root  comer  frequency  scaling.  However,  most  of  the  events  considered  have  local  magnitudes 
less  than  3.0,  so  the  comer  frequencies  are  near  or  beyond  the  upper  end  of  the  frequency  band 
inverted.  The  results  are  therefore  not  very  sensitive  to  our  source  parameterization. 

All  of  the  results  discussed  previously  were  obtained  assuming  the  simplified  Mueller-Murphy 
(1971)  explosion  source  function  (3.1.4)  with  cube-root  corner  frequency  scaling  (3.1.6).  We 
also  inverted  the  Lg  data  assuming  the  Brune  (1970,  1971)  earthquake  source  model  (3.1.5) 
with  cube-root  corner  frequency  scaling.  For  the  same  Qq,  the  Brune  source  model  leads  to  Q 
with  stronger  frequency  dependence.  For  example,  when  Lg  Qq  is  350,  T)  is  0.41  for  the 
Mueller-Murphy  model  and  0.48  for  the  Brune  model.  The  reason  for  this  slight  difference  is 
that  the  Brune  model  has  a  less  abrupt  corner  frequency  (Figure  3.1),  and  the  reduced  high  fre¬ 
quency  source  contribution  is  compensated  for  by  a  higher  Q  at  high  frequency.  The  Brune 


model  inversions  give  a  data  variance  that  is  2%  greater  than  that  obtained  with  the  Mueller- 
Murphy  model. 

If  we  assume  a  flat  source  (no  comer  frequency),  the  results  are  nearly  the  same  as  those  hum 
the  Mueller-Murphy  (1971)  model,  with  a  slight  increase  in  the  data  variance.  This  mesely 
indicates  that  the  comer  frequencies  of  most  of  the  events  are  greater  than  7  Hz.  We  have  nut 
explored  other  alternative  source  models,  but  this  result  indicates  that  ourLg  attenuation  results 
are  not  sensitive  to  the  details  of  the  source  parameterization. 

The  Pn  inversion  is  more  sensitive  to  source  parameterization  than  the  Lg  inversion  as  a  result 
of  both  radiation  pattern  and  the  increased  bandwidth.  We  used  the  Mueller-Murphy  (1971) 
model  for  all  Pn  inversions.  We  have  not  run  the  Pn  inversion  with  the  Brune  (1970,  1971) 
earthquake  source  function,  although  we  would  expect  a  response  similar  to  that  observed  for 
Lg.  That  is,  a  slight  increase  in  the  frequency  dependence  of  Q  to  compensate  for  the  reduced 
high  frequency  source  contribution.  We  have  not  analyzed  the  sensitivity  of  our  estimate  with 
respect  to  high  frequency  source  falloff,  but  it  is  clear  that  an  increase  in  Q  must  accompany 
an  increased  source  spectral  decay  to  preserve  the  fit  at  high  frequency.  Our  results  for  Pn 
must  therefore  be  considered  relative  to  an  (D2  source. 

6  3.2  Geometric  spreading 

The  geometric  spreading  rate  for  Lg  is  well-constrained  both  empirically  and  theoretically,  and 
our  assumptions  are  described  in  Section  3.1.  All  of  the  data  used  in  this  study  are  at  ranges 
greater  than  100  km.  Therefore,  r0  (in  3.1.2)  trades  off  directly  with  inverted  moment  and  has 
no  effect  on  the  Q^gif)  estimate. 

The  geometric  spreading  rate  of  Pn  is  much  more  difficult  to  estimate  with  confidence.  Results 
presented  earlier  were  obtained  by  assuming  r-13  Pn  spreading.  In  this  section,  we  examine 
the  results  for  alternatives,  in  particular  r~l  (spherical  spreading)  and  r_1  5  were  used.  Table 
6.3  describes  the  resulting  Pn  models,  all  required  to  give  a  range-independent  So{Lg)lS0{Pn) 
ratio  for  the  Lg  Q  model  of  Section  5  (#4  in  Table  6.1).  As  the  Pn  spreading  rate  increases, 
Qo(Pn)  increases,  T)  decreases,  and  Mq  increases.  The  range-independent  value  of  the  long 
period  source  level  ratio  can  be  used  to  extract  the  relative  Lg  excitation  of  explosions  and 
earthquakes  (5.3.1).  The  K  in  Table  6.3  indicates  the  relative  excitation  of  Lg  for  explosions 
and  earthquakes.  Note  that  assuming  r'1  Pn  spreading  results  in  greater  Lg  excitation  for 
explosions  than  for  earthquakes  of  equal  moment  (tc  >  1).  This  contradicts  empirical 


Table  6.3.  Pn  parameter  trade-off  as  a  function  of  spreading  rate 
Spreading  Rate  Q0  r)  c  LogM **p  k 


246 

0.53 

18.6 

1.09Ml  +  16.95 

1.11 

0.0623 

300 

0.49 

29.2 

1.12 Ml  +  17.57 

0.26 

0.0637 

352 

0.45 

68.3 

1.14M/  +  17.99 

0.09 

0.0651 

observations  and  so  argues  for  a  Pn  spreading  rate  greater  than  r-1  On  the  other  hand,  a  Pn 
spreading  rate  of  r-1  5  gives  K  =  0.09,  or  11  times  greater  Lg  for  earthquakes  than  explosions. 
This  is  much  too  large  a  difference.  The  earthquake  Lg  excitation  for  r"1  3  Pn  spreading  is  4 
times  the  explosion  excitation,  consistent  with  the  observations  of  Willis  (1963)  and  Pomeroy 
(1977).  This  is  also  a  result  consistent  with  theoretical  simulations  of  Pn  propagation  in  realis¬ 
tic  earth  models  (Wallace,  personal  communication). 

6.4  Systematic  variations  in  source  moment  versus  Lg  amplitude 

We  have  not  accounted  for  any  azimuthal  variations  in  Q,  but  have  combined  data  from  all 
azimuths  into  an  inversion  for  a  single,  average  Q  model.  Examination  of  the  fit  of  the  model 
to  the  data  (Appendix  B)  demonstrates  the  validity  of  this  procedure.  However  we  do  note 
some  systematic  variations  in  the  Lg  amplitude  (represented  by  M i)  for  fixed  moment  obtained 
from  the  inversion.  In  particular,  in  the  Mq-M ^  data  in  Figure  5.5  there  is  a  set  of  explosions 
with  lower  magnitudes  for  fixed  moment  than  is  indicated  by  the  general  trend  of  the  data. 
These  are  explosions  from  mines  K1-K12  at  azimuths  between  34°  and  46°  and  from  mines 
V1-V13  at  azimuths  between  57°  and  84°  (Figure  4.3).  The  neighboring  mines,  E1-E9 
between  91°  and  97°,  are,  however,  consistent  with  the  general  trend.  Figure  6.4  displays  the 
Lg  explosion  moments  versus  magnitude  for  the  V  mines,  K  mines,  and  E  mines  along  with 
the  curve  described  by  (5.4.2).  In  general,  events  from  the  V  and  K  mines  have  smaller  Lg 
amplitudes  for  a  given  moment  than  events  from  E  mines.  Figure  6.5  compares  NORESS  verti¬ 
cal  component  recordings  of  an  E7  mine  blast  to  recordings  of  a  V8  mine  blast.  Both  events 
have  NORESS  magnitudes  of  3.3.  The  disparate  PnlLg  ratios  for  the  two  events  indicates  that 
either  Lg  propagation  is  structurally  inhibited  north  of  E7  or  that  different  mining  practices 
result  in  distinctly  different  Pn  to  Lg  excitation  ratios.  If  Lg  propagation  is  structurally  inhi¬ 
bited,  in  must  be  a  broad  band  blockage  because  spectra  from  V8  and  E7  mines  are  consistent 
with  the  same  Q  model  (Appendix  B). 

6.5  Anomalous  events 

Appendix  A  lists  events  that  were  processed  but  not  included  in  the  inversion.  Most  of  those 
were  excluded  because  of  low  signal  to  noise  ratios,  signals  from  multiple  events  were  mixed, 
or  they  had  uncertain  locations  due  to  ambiguous  phase  association.  Also,  there  were  1 1 
events  that  occurred  within  a  24  hour  time  period  off  the  southwest  coast  of  Norway  that  were 
not  included  because  their  spectra  were  inconsistent  with  those  from  events  at  similar  distances. 
The  event  epicenters  are  plotted  in  Figure  A.l  (events  15-23  and  25-26  of  Table  A.I).  We  are 
not  certain  whether  these  are  off-shore  earthquakes  or  underwater  explosions  (Henson  and 
Bache,  1986).  Figure  6.6  is  a  comparison  of  Lg  spectra  for  a  Titania  mine  blast  and  of  one  of 
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Figure  6.4.  Lg  inverted  explosion  moment  versus  NORESS  local  magnitude  for  model  4  of 
Table  6.1.  Only  explosions  from  K,  V,  and  E  mines  (Table  4.2)  are  included. 
The  triangles  are  K-mine  explosions,  asterisks  are  V-mine  explosions,  and  cir¬ 
cles  are  E-mine  explosions.  The  line  is  the  LogA/g-magnitude  relation  derived 
using  all  of  the  explosions  in  the  data  set  (S.4.2). 
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Figure  6.6.  Comparison  of  Lg  spectra  of  a  Titania  mine  blast  (event  63,  Table  4.1)  and  an 
off-shore  event  (event  20,  Table  A.l).  The  off-shore  event  is  at  approximately 
the  same  azimuth  and  only  80  km  further  from  NORESS  than  the  Titania  mine. 
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the  off-shore  events.  The  increased  high  frequency  decay  of  the  off-shore  event  spectrum  can¬ 
not  be  accounted  for  by  the  different  path  lengths.  The  Pn  spectra  of  the  two  events  have  a 
similar  relationship.  There  are  several  possible  explanations  for  the  difference  between  the 
spectra.  If  the  off-shore  events  are  explosions,  their  high  frequency  energy  could  be  depleted 
by  near-source  sediment  resonance.  One  might  also  speculate  that  the  difference  is  caused  by 
some  structural  feature  between  the  off-shore  events  and  the  coast  of  Norway.  However,  we 
note  that  events  66  and  97  of  Table  4.1  are  off-shore  southwest  Norway  events  that  have  spec¬ 
tra  that  are  consistent  with  the  Titania  mine  blast  spectra  (Appendix  B).  This  suggests  that  the 
difference  is  either  from  source  differences  or  near-source  structure. 


7.  SIMULATION  OF  REGIONAL  Pn  AND  Lg  SPECTRA 


We  have  developed  an  accurate  parameterization  of  the  observed  spectra  of  regional  events 
recorded  at  N'ORESS  in  terms  of  familiar  seismological  parameters  characterizing  the  source 
and  range  dependence  of  these  spectra.  Using  these  parameters  we  can  "predict"  the  spectra  for 
a  chosen  source  and  range  and  be  confident  that  it  is  correct  with  well-defined  uncertainty  of 
our  parameters.  Our  confidence,  of  course,  degrades  if  we  choose  a  source  or  range  outside 
the  bounds  of  our  experience;  that  is,  if  we  attempt  to  extrapolate  our  results  to  situations  not 
yet  encountered. 

Figure  7.1  displays  predicted  range-dependent  Pn  spectra,  using  the  inversion  results  of  Section 
5,  for  a  magnitude  3.0  explosion.  The  relevant  parameters  are  Q(f)  =  300/ 049,  M0  = 
8.6  x  10“°  dyne-cm  (50  =  21.9),  fc  =  10.4  Hz,  and  a  Pn  spreading  rate  of  r~1,3.  Note  that 
these  results  do  not  depend  on  resolving  the  parameter  trade-offs;  all  of  the  models  identified 
in  Section  6  will  give  similar  predictions.  Also  shown  is  the  average  N'ORESS  noise  estimated 
from  samples  taken  prior  to  Pn  for  many  events  (Henson  and  Bache,  1986).  The  Pn  spectra 
approximately  parallel  the  noise  curve  at  a  distance  of  400-500  km,  converging  with  it  at 
longer  ranges  and  diverging  from  it  at  shorter  ranges.  This  observation  is  supported  by  average 
spectral  density  curves  constructed  by  Ringdal,  et  al.  (1986)  using  the  high  frequency  element 
of  the  NORESS  array.  Figure  7.2  is  the  corresponding  plot  for  Lg.  The  Lg  parameters  are 
Q(J)  =  350/041,  Mq  =  9.8  x  1020  dyne-cm  (Sq  =  20.7),  fc  =  10.4  Hz,  and  cylindrical  spread¬ 
ing  with  r0  =  100  km.  The  Lg  spectrum  reaches  the  noise  level  at  a  much  lower  frequency 
than  the  Pn  spectrum,  while  at  long  periods  the  Lg  amplitude  exceeds  the  Pn  amplitude  by  as 
much  as  a  factor  of  10.  This  indicates  that  these  regional  seismograms  are  characterized  by  Lg 
the  largest  amplitude  phase  and  Pn  having  a  higher  dominant  frequency. 

The  largest  uncertainty  in  the  "predicted"  spectra  in  Figures  7.1  and  7.2  is  the  absolute  level 
for  the  magnitude  chosen;  that  is,  the  scatter  in  the  moment-magnitude  relation  (Figure  5.5). 
We  must  also  account  for  the  standard  error  in  the  spectral  estimate  for  a  particular  moment 
(Figure  6.1),  but  this  is  relatively  small.  Combining  variances  for  these  two  contributions  to 
the  uncertainty,  we  compute  the  standard  error  for  our  spectral  estimate.  This  is  illustrated  in 
Figure  7.3  with  the  predicted  Pn  spectrum  at  800  km  from  a  NORESS  magnitude  3.0  explo¬ 
sion.  The  Lg  spectral  estimate  has  lower  variance  for  fixed  moment  (Figure  6.1),  but  the  vari¬ 
ance  of  the  Lg  explosion  moment-magnitude  relation  is  larger  than  the  corresponding  variance 
for  Pn.  The  result  is  that  the  standard  error  in  our  predicted  Lg  spectra  is  about  the  same  as 
the  error  for  Pn,  which  is  illustrated  in  Figure  7.3. 


Predicted  NORESS  Pn  displacement  spectra  at  3  epicentral  distances  for  a  mag¬ 
nitude  3.0  explosion,  based  on  the  inversion  results.  The  average  NORESS  noise 
was  estimated  from  samples  taken  prior  to  Pn  for  many  events  (Henson  and 
Bache,  1986). 
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Figure  7.2.  Predicted  NORESS  Lg  displacement  spectra  at  3  epi central  distances  for  a  mag¬ 
nitude  3.0  explosion,  based  on  the  inversion  results.  The  average  NORESS  noise 
was  estimated  from  samples  taken  prior  to  Pn  for  many  events  (Henson  and 
Bache,  1986). 


8.  CONCLUSIONS 


This  report  summarizes  results  of  a  detailed  analysis  of  regional  wave  attenuation  along  con¬ 
tinental  paths  to  the  NORESS  array  in  Norway.  We  have  developed  and  implemented  a  gen¬ 
eralized  inversion  of  log  amplitude  spectra  that  simultaneously  estimates  seismic  moment  and 
apparent  attenuation.  The  method  uses  both  the  spatial  and  spectral  decay  of  observed  signal 
amplitudes  to  separate  source  and  path  contributions.  The  inversion  was  applied  to  Pn  and  Lg 
spectra  from  190  regional  events  recorded  at  NORESS  over  distances  of  200  to  1300  km. 
Based  on  adequate  signal/noise,  the  Lg  spectra  were  inverted  between  1  and  7  Hz  and  Pn 
spectra  between  1  and  15  Hz.  The  main  conclusions  of  the  study  are  as  follows: 

•  Data  from  all  190  regional  events  are  adequately  modeled  by  a  simple  Or  source 
and  a  single  frequency -dependent  Q  model.  Theoretical  spectra  derived  from  this 
simple  parameterization  were  compared  to  322  observed  regional  wave  spectra.  The 
ability  of  the  model  to  reproduce  the  important  spectral  characteristics  of  such  a 
large  number  and  variety  of  observed  data  provides  support  for  the  results  that  goes 
well  beyond  that  in  most  published  studies  of  regional  wave  attenuation. 

•  Regional  Lg  spectra  are  adequately  represented  from  1  to  7  Hz  by  a  power-law  fre¬ 
quency  dependence  of  apparent  attenuation  given  by  Q(f)  -  350/ 041.  Our  short 
time  windows  include  the  onset  of  Lg  which  primarily  consists  of  higher  mode  sur¬ 
face  waves  sampling  the  lower  crust.  Our  Qigif)  is  therefore  an  approximation  to 
the  Q  of  the  lower  crust,  although  it  includes  the  effects  of  apparent  attenuation  due 
to  scattering. 

•  Amplitudes  of  regional  Pn  phases  are  more  sensitive  to  source  radiation  pattern, 
focussing  and  defocussing,  and  scattering.  Nevertheless,  a  successful  parameteriza¬ 
tion  of  regional  P  wave  attenuation,  particularly  at  high  frequency,  has  important 
applications  in  treaty  monitoring  seismology.  We  have  found  that  regional  Pn  spec¬ 
tra  of  152  events  are  consistent  with  r_1  3  geometric  spreading  and 
Qpnif)  =  300/  049.  We  do  not  attempt  to  distinguish  intrinsic  absorption  from 
scattering  and  explicitly  acknowledge  that  our  Q(f)  is  an  empirical  parameterization 
of  the  data. 


•  Our  estimated  seismic  moments  as  a  function  of  local  magnitude  are  generally  con¬ 
sistent  with  near-fteld  studies.  In  that  magnitude  does  not  enter  the  inversion,  this 
lends  considerable  support  to  the  derived  Q  model.  Comer  frequencies,  while  in 
agreement  with  near-field  studies,  are  not  clearly  resolved  by  our  data. 

•  For  a  fixed  set  of  source  and  spreading  assumptions,  the  inversion  defines  a  broad 
minimum  in  the  data  residuals  corresponding  to  a  suite  of  models  that  fit  the  data 
equally  well  in  a  least-squares  sense.  However,  an  important  constraint  on  the 
attenuation  estimates  is  the  consistency  of  the  derived  source  parameters  from  the 
independently  inverted  Pn  and  Lg  spectra.  This  provides  a  firm  internal  consistency 
check  on  our  attenuation  results  and  reduces  the  trade-offs  among  model  parameters 
to  a  set  of  model  pairs.  That  is,  for  fixed  Qigif),  the  source  constraint  places  tight 
bounds  on  Qpn(f).  A  trade-off  between  Mq-M l  variance  and  data  variance  is  used 
to  define  the  range  of  acceptable  model  pairs.  The  "preferred  model"  was  merely 
selected  from  the  mid-range  of  these  acceptable  solutions.  The  range  of  acceptable 
Q  models  for  Lg  is  Qq  between  300  and  400  and  r|  between  0.46  and  0.37.  Simi¬ 
larly,  for  Pn  with  r_1  3  spreading,  the  range  is  Q0  between  250  and  375  with  T| 
between  0.53  and  0.44. 

•  The  Qigif)  results  are  insensitive  to  the  details  of  our  source  parameterization.  Most 
events  considered  had  local  magnitudes  less  than  3.0,  and  corner  frequencies  greater 
than  the  upper  limit  of  the  bandwidth  inverted.  The  sensitivity  of  Qpn(f)  to  source 
parameterization  is  unknown,  although  a  spectral  decay  faster  than  or  must  be 
accompanied  a  higher  QPn  to  preserve  the  fit  at  high  frequency.  Therefore,  the 
derived  Qpn  function  must  be  considered  relative  to  an  to2  source  model. 

•  Pn  geometric  spreading  rates  of  r~\  r~13,  and  r-1  5  were  investigated.  In  each 
case,  the  inversion  resulted  in  a  Q  model  that  could  adequately  reproduce  the  Pn 
spectra.  However,  only  r-1  3  resulted  in  an  earthquake  to  explosion  Lg  excitation 
ratio  which  is  consistent  with  other  empirical  studies. 

•  Data  from  all  azimuths  were  combined  into  an  inversion  for  a  single,  average  Q 
model.  No  evidence  for  an  azimuthal  dependence  in  spectral  decay  was  discovered, 
although  explosions  from  different  mines  were  found  to  produce  distinctly  different 
observed  Lg  to  Pn  ratios.  It  is  not  clear  whether  the  low  Lg/Pn  ratios  are  due  to 
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structural  inhibition  of  Lg  propagation  or  to  different  mining  practices  resulting  in 
different  Pn  to  Lg  excitation  ratios. 

•  The  Lg  and  Pn  attenuation  models  are  used  to  predict  range-dependent  spectral 
amplitudes  for  events  of  arbitrary  magnitude.  For  a  broad  distance  range,  the  model 
predicts  regional  seismograms  characterized  by  Lg  the  largest  amplitude  phase  and 
Pn  having  a  higher  dominant  frequency. 

The  final  result  of  the  inversion  is  an  accurate  parameterization  of  observed  amplitude  spectra 
of  regional  events  recorded  at  NORESS  that  can  be  used  to  address  a  number  of  seismological 
problems  related  to  wave  propagation  in  the  region  and  to  the  treaty  monitoring  capabilities  of 
small  regional  networks. 
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APPENDIX  A.  PROCESSED  EVENTS  NOT  INCLUDED  IN  THE  INVERSION 


Section  4  tabulates  the  events  used  in  the  generalized  inversion.  In  addition  to  these,  there  were 
68  events  processed  that  were  not  included  in  the  inversion.  Table  A.l  identifies  this  set  of 
events.  Most  of  these  events  were  excluded  because  of  low  signalnoise  ratios  Others  were 
excluded  because  they  were  out  of  the  regional  distance  range  to  NORESS.  Event  6  of  Table 
A.l  is  a  regional  Soviet  Peaceful  Nuclear  Explosion  (PNE)  with  mb  =  5.0.  This  event  was 
excluded  because  it  was  considerably  larger  than  other  events  in  the  data  set  and  because  its 
source  spectrum  is  likely  to  be  inherently  different  from  spectra  of  the  small,  near-surface, 
mine  blasts  comprising  the  majority  of  the  events  used. 

Many  of  the  events  in  Table  A.l  were  processed  to  test  and  evaluate  the  performance  of  die 
"expert  system"  currently  under  development  at  SAIC.  These  events  were  selected  because 
they  are,  in  some  sense,  difficult  as  a  result  of  low  signal  noise  ratios,  mixed  events,  problems 
in  location  due  to  phase  misassociation  ,  etc.  A  few  of  the  events  listed  in  Table  A.l  were 
excluded  because  signals  from  multiple  events  were  mixed.  Others,  without  independent  net¬ 
work  solutions,  had  ambiguous  phase  associations  and  therefore  uncertain  locations.  There 
were  a  couple  of  events  that  were  reported  by  an  independent  network,  but  that  location  did 
not  agree  with  the  arrival  pattern  recorded  at  NORESS. 

More  fundamentally,  however,  there  were  16  events  that  had  regional  wave  spectra  inconsistent 
with  the  other  events  at  approximately  the  same  distance  These  are  events  15-23,  25-26,  40, 
45,  and  50-52  of  Table  A.l.  Their  locations  are  plotted  on  a  map  relative  to  NORESS  in  Fig¬ 
ure  A.l.  For  each  of  these  events,  the  Pn  and  Lg  specrra  decay  with  frequence  considerably 
faster  than  is  predicted  by  the  model  derived  from  the  events  in  Table  4.1.  Event  50  is  located 
only  3  km  from  mine  NYG  (Table  4.2)  and  the  other  15  events  are  all  presumed  euithquakes. 
Event  45  does  not  have  an  independent  network  solution,  although  there  is  nothing  apparent!) 
wrong  with  the  ROKAPP  phase  association.  Events  40,  45.  and  50  appear  to  he  anomalous 
events"  in  that  other  events  along  the  same  azimuth  and  approximate  distance  have  spectra 
consistent  with  the  majority  of  the  events  studied.  Regional  spectra  from  events  51  and  52 
(located  near  the  northern  tip  of  Denmark)  may  reflect  greater  path  attenuation  however  the 
interpretation  is  ambiguous  Events  15-23  and  25-26  all  occurred  within  a  24  hour  pcnod  and 
were  located  off  the  southwest  coast  of  Norway  (Figure  A.l)  It  is  noi  known  whether  these 
are  off-shore  earthquakes  or  underwater  explosions  (Henson  and  Ruche.  lOXhi  Because  (ties 
are  so  nearly  <onfemporuneous  and  of  similar  magnitude,  it  is  difficult  to  rule  :  anomalous 
source  characteristics  as  contributing  to  the  difference  between  these  s|\\ii.i  ,u  1  ih<-ce  of  the 


Table  A.l.  Events  processed  not  included  in  the  inversion. 
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5- 17-85 

6- 18-85 

6- 30-85 

7- 11-85 
7-18-85 
7-25-85 

10- 30-85 

11-  6-85 
11-  8-85 
11-14-85 
11-14-85 
11-15-85 
11-15-85 
1 1-20-85 
11-20-85 
11-20-85 
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11-20-85 
11-20-85 
11-20-85 
11-21-85 
11-21-85 
11-21-85 
11-21-85 
11-21-85 

11- 29-85 

12-  4-85 
12-11-85 
12-13-85 
12-13-85 
12-16-85 
12-13-85 
12-23-85 
12-26-85 

I-  3-86 
1-  9-86 
1-  9-86 
1-13-86 
1-19-86 
1-23-86 
1-24-86 

1- 31-86 

2-  3-86 
2-13-86 


0:57  49.90N  79.00E 

0:25  55.00N  93.00E 

3:58  55.00N  91.00E 

2:39  49.90N  79.00E 


21:15 

3:11 

7:40 

9:00 

0:08 

10:44 

18:20 

11:00 
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22:11 
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22:57 

23:11 
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14:07 

14:18 

14:48 

15:05 
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12:24 

14.00 

20:21 
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14:14 

16:44 

2:35 
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13:33 

6:16 

9:18 

20:35 

9:36 

4:59 


66.00N 

49.90N 

72.08N 

64.70N 

51.68N 

62.70N 

59.40N 

67.60N 

63.20N 

57.61N 

57.66N 

57.63N 

57.66N 

57.69N 

57.64N 

57.58N 

57.75N 

57.06N 

54.80N 

57.44N 

58.30N 

59.60N 

60.73N 

15.57N 

64.70N 

67.60N 

67.  ION 

60.38N 

50.21N 

58.90N 

59.56N 

54.70N 

51.25N 

64.60N  : 

64.95N 

50.21N 

64.70N 

64.70N 

58.46N 

60.30N 


40.90E 
78.20E 
1.34W 
30.70E 
16.20E 
17.76E 
34.20E* 
34.20E 
27.80E 
5.67E 
:  5.72E 
6.27E 
5.35E 
5.45E 
5.62E 
5.49E 
5.30E 
6.36E 
6.50E* 
5.74E 
4.80E* 
18.70E* 
11.31E 
16.16E 
30.70E 
34.00E 
20.60E 
1.90E 
12.39E 
33.40E* 
7.45E 
19.50E* 
15.57E 
20.40E* 
12.13E 
I2.40E 
30  70E 
30.70E 
12  10E 
15.00E* 


EX 

EX-PNE 

EX 

EQ 

EX-V10 

EQ 

EQ 

EQ 

EX-K5 

EX-M7 

EQ 

EQ 

EQ 

EQ 

HQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EX-V10 

EX-K2 

EX-R1 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EQ 

EX-V10 

EX-V10 

EQ 

EQ 


TO® 


Event 

Date 

Time 

Location 

Type 

Magnitude 

46 

2-16-86 

15:04 

58.29N  30.40E 

EQ 

<2.0 

(H) 

47 

2-24-86 

1:13 

59.30N  3.70E 

EQ 

1.5 

(R) 

48 

4-  8-86 

10:06 

59.30N  27.60E 

EX-E6 

2.1 

(H) 

49 

4-  8-86 

10:30 

64.33N  20.56E* 

EQ 

2.4 

(H) 

50 

4-  8-86 

1 1:28 

60.39N  5.29E 

— 

1.9 

(B) 

51 

4-  9-86 

8:16 

57.31N  10.60E 

EQ 

2.5 

(B) 

52 

4-  9-86 

8:27 

57.30N  10.48E 

EQ 

2.5 

(B) 

53 

4-18-86 

0:44 

59.22N  1.42E 

EQ 

2.4 

(B) 

54 

4-30-86 

6:05 

57.28N  6.19E 

EQ 

2.2 

(B) 

55 

4-30-86 

6:23 

57.40N  6.17E 

EQ 

2.3 

(B) 

56 

5-  2-86 

10:45 

60.35N  5.08E 

— 

1.8 

(B) 

57 

5-  5-86 

15:37 

60.24N  5.28E 

— 

1.5 

(B) 

58 

5-27-86 

18:36 

6 1.5 IN  5.84E 

EQ 

1.6 

(B) 

59 

6-11-86 

14:08 

59.39N  18.66E 

EQ 

2.1 

(H) 

60 

6-17-86 

12:03 

58.59N  18.14E 

EQ 

2.6 

(H) 

61 

7-19-86 

11:54 

54. 10N  21.67E 

EQ 

2.3 

(H) 

62 

7-30-86 

13:43 

59.50N  30.00E 

EX 

2.0 

(H) 

63 

7-30-86 

22:30 

59.77N  11.05E 

EQ 

2.4 

(B) 

64 

8-10-86 

5:01 

59.99N  5.34E 

— 

1.5 

(B) 

65 

8-18-86 

8:25 

59.50N  25.00E 

EX-E3 

2.0 

(H) 

66 

10-12-86 

14:53 

60.92N  5.1  IE 

EQ 

1.9 

(B) 

67 

10-30-86 

9:21 

59.60N  10.80E* 

EQ 

1.2 

(R) 

68 

10-30-86 

9:22 

60.80N  6.40E* 

EQ 

1.6 

(R) 

*  RONAPP  Location 
t  Alewine,  1985 
(B)  Bergen  network  magnitude 
(H)  Helsinki  network  magnitude 
(N)  NORSAR  magnitude 
(P)  PDE  magnitude  (Mh) 

(R)  RONAPP  magnitude 


Events  With  Anomalous  Spectra 


Figure  A.l  Locations  of  events  with  anomalous  spectra.  The  events  are  identified  by  the 
event  numbers  in  Table  A.l.  Events  15-23  and  25-56  all  occurred  within  a  24 
hour  time  period  and  may  be  underwater  explosions  (Henson  and  Bache,  1986). 


nearby  Titania  mine.  Clearly,  the  other  possibility  is  considerable  path  attenuation  between 

these  event  epicenters  and  the  coast  of  Norway.  These  events  will  be  discussed  in  more  detail 
in  Section  6.5. 


APPENDIX  B.  INVERSION  RESULTS 


Regional  Pn  and  Lg  spectra  of  190  events  (Table  4.1)  were  inverted  for  Q(f),  seismic  moment, 
and  the  parameter,  c,  relating  corner  frequency  to  long  period  source  level.  The  source  was 
parameterized  by  the  Mueller-Murphy  explosion  model  (3.1.4)  and  geometric  spreading  was 
modeled  by  (3.1.2)  with  m  =  1/2  and  r0  =  100  km  for  Lg  and  m  =  1.3  with  r0  =  1  km  for  Pn. 
In  terms  of  a  power-law  frequency  dependence,  the  path  result  for  Pn  is  Q(f)  =  300/  0  jg  and 
for  Lg  is  Q(f)  =  350/°  41 .  The  parameter,  c,  derived  from  Pn  spectra  is  29.2  and  from  Lg 
spectra  is  28.7. 


Table  B.l  lists  the  results  of  the  inversion  for  source  parameters.  The  event  numbers 
correspond  to  those  of  Table  4.1  and  the  mine  identifications  for  the  explosions  are  given  in 
Table  4.2.  Comer  frequencies  were  not  parameters  of  the  inversion  for  each  event,  but  were 
calculated  from  (3.1.6)  using  the  inversion  results  for  c.  For  the  explosions,  seismic  moment 
was  estimated  from  Sq(Ph)  using  (3.1.7)  assuming  a  surface  compressional  velocity  of  5  km/s 
and  a  density  of  2.5  gm/cm3.  The  earthquake  moments  were  estimated  from  S0(Lg )  using 
(3.1.9)  with  an  average  crustal  shear  wave  velocity  of  3.5  km/s  and  a  crustal  density  of  2.7 
gm/cm  .  Because  the  radiation  patterns  are  unknown,  moments  were  not  estimated  for  earth¬ 
quakes  for  which  only  Pn  spectra  were  included  in  the  inversion  or  for  events  of  unknown 
source  type.  Similarly,  moments  were  not  estimated  for  explosions  without  Pn  spectra  because 
of  the  uncertainty  in  the  relative  Lg  to  Pn  source  excitation  (3.1.10). 

Figure  B.l  displays  theoretical  and  data  spectra  for  all  of  the  events  used  in  the  inversion.  The 
theoretical  spectra  were  computed  using  (3.1.1)  with  the  results  of  the  inversion  for  Q(f)  and 
S(J).  Each  spectrum  is  identified  by  the  event  number  in  Table  4.1.  The  spectra  have  been 
vertically  offset  for  display  purposes, 
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Table  B.l  Inverted  source  parameters 


92 

3.1 

EX-E4 

.... 

93 

2.6 

EX-V10 

26.84 

94 

3.1 

EX-K2 

57.25 
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Figure  B.l 


Theoretical  Pn  and  Lg  spectra  compared  to  observed  spectra  for  the  190  events 
used  in  the  generalized  inversion.  The  theoretical  spectra  were  computed  using 
the  inversion  parameters  described  in  Appendix  B  and  listed  in  Table  B.l. 
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Figure  B.l  (Continued). 
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